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Global  prediction  of  the  ocean  weather  (e.g.,  surface  mixed  layer,  meandering 
currents  and  fronts,  eddies,  and  eoastally  trapped  waves)  has  been  feasible  onl\ 
sinee  the  turn  of  the  eentury  w  hen  sufficient  computing  pow  er  and  real-time  data 
became  available.  Satellite  altimetry  is  the  key  observing  system  for  mapping  ocean 
eddies  and  current  meanders,  but  sea  surfaee  temperature,  temperature  and  salinity 
profiles,  and  atmospheric  forcing  arc  also  essential.  The  multinational  Global  Ocean 
Data  Assimilation  Experiment  (GODAE)  has  fostered  the  development  of  eddy- 
permitting  and  eddy-resolving  basin-scale  to  global  prediction  systems  in  several 
countries.  In  this  work,  the  foeus  is  on  eddy-resolving  global  ocean  prediction,  the 
value  of  ocean  model  skill  in  dynamical  interpolation  of  data  and  in  forecasting,  plus 
the  increased  skill  gained  from  using  eddy-resolving  versus  eddy-permitting  models. 
In  ocean  prediction,  one  must  consider  the  classes  of  ocean  response  to  atmospheric 
forcing  for  different  phenomena  and  regions.  In  oeean  noweasting,  this  is  important 
in  considering  the  relative  impaet  of  the  atmospheric  forcing  and  different  oeean 
data  types.  In  forecasting,  these  classes  impact  whether  or  not  the  timeseale  for 
oceanie  predietive  skill  is  limited  by  the  timescale  for  atmospheric  predictive  skill 
or  whether  a  coupled  oeean-atmosphere  model  would  be  advantageous.  Results 
from  the  existing  eddy-resolving  global  oeean  prediction  systems  demonstrate 
forecast  skill  up  to  1  month,  globally  and  over  most  subregions.  Outside  surface 
boundary  layers  and  shallow  water  regions,  the  forecast  skill  typically  is  only 
modestly  impacted  by  reverting  toward  elimatological  forcing  after  the  end  of  the 
atmospheric  forecast  versus  using  analysis-quality  forcing  for  the  duration. 


I  INTRODUCTION 


Ocean  Modeling  in  an  Eddying  Regime 
Geophysical  Monograph  Series  1 77 
Copyright  2008  by  the  American  Geophysical  l  nion 
10. 1029/1 77GM  21 


Adequate  real-time  data  input,  computing  power,  numeri¬ 
cal  oeean  models,  data  assimilation  capabilities,  atmospheric 
forcing,  and  bottom  topography  are  key  elements  required 
for  sueeessful  eddy-resolving  global  prediction  of  the  “ocean 
weather",  such  as  the  surfaee  mixed  layer,  equatorial  and 
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coastally  trapped  waves,  upwelling  of  cold  water,  meander¬ 
ing  of  currents  and  fronts,  eddies,  Rossby  waves,  and  the 
associated  temperature,  salinity,  currents,  and  sea  surface 
height  (SSH).  Only  since  the  turn  of  the  century  have  all 
the  key  elements  finally  reached  a  status  that  makes  global 
ocean  weather  prediction  feasible.  The  focus  of  this  chapter 
is  on  the  roles  that  eddy-resolving  ocean  models  perform  in 
accurate  nowcasting  (estimating  the  present  state)  and  fore¬ 
casting  of  the  ocean  weather  on  timescales  up  to  a  month. 
However,  much  of  the  content  is  also  relevant  to  climate 
prediction. 

In  the  remainder  of  the  introduction  the  relevance  of  this 
chapter  to  climate  is  briefly  discussed  in  section  1.1.  Section  1 .2 
outlines  the  history  of  the  first  1 0  years  of  basin-scale  to  glo¬ 
bal  ocean  weather  prediction,  the  product  of  a  multinational 
effort  largely  fostered  by  the  Global  Ocean  Data  Assimila¬ 
tion  Experiment  (GODAE).  In  section  1.3,  ocean  prediction 
is  discussed  in  relation  to  four  classes  of  ocean  response  to 
atmospheric  forcing,  how  they  impact  ocean  prediction,  and 
how  they  affect  the  involvement  of  ocean  models.  Section  2 
presents  results  illustrating  the  roles  that  ocean  models  play 
in  nowcasting  and  forecasting  with  emphases  on  results  from 
eddy-resolving  global  ocean  models  and  prediction  systems 
and  the  relation  to  classes  of  responses  to  atmospheric  forc¬ 
ing  These  include  dynamical  interpolation  of  assimilated 
data,  the  value  of  atmospheric  forcing  as  a  data  ty  pe,  sensi¬ 
tivity  to  ocean  model  resolution  (e.g.,  eddy-resolving  versus 
eddy-permitting  and  impact  on  the  space  scales  that  can  be 
mapped  using  nadir  beam  satellite  altimetry),  use  of  mean 
SSH  from  an  eddy-resolving  global  ocean  model  in  creating 
a  mean  SSH  field  for  addition  to  the  anomalies  from  satellite 
altimetry,  and  the  use  of  simulated  data  from  eddy-resolving 
ocean  models  in  testing,  evaluating,  and  improving  data  as¬ 
similation  skill.  The  latter  is  also  an  example  of  observing 
system  simulation  and  assessment. 

1 ,  /.  Relevance  to  Climate 

To  date,  multidecadal  climate  prediction  models  and  sea¬ 
sonal  to  interannual  climate  forecasts  (focused  primarily  on 
El  Nino/La  Nina  conditions)  have  used  coarse-resolution 
ocean  models  that  are  unable  to  depict  realistic  ocean  fronts. 
As  climate  predictions  become  increasingly  region-specific, 
eddy-resolving  ocean  models  will  be  required  for  success  in 
many  regions,  for  example,  for  accurate  positioning,  strength 
and  eastward  penetration  of  major  ocean  currents  and  fronts. 
As  specific  examples,  eddy-resolving  ocean  models  are  re¬ 
quired  for  realistic  prediction  of  the  Kuroshio/Oyashio  cur¬ 
rent  system  in  the  western  North  Pacific  and  the  Gulf  Stream 
current  system  in  the  western  North  Atlantic.  In  the  Pacific, 
such  models  must  simulate  a  strongly  inertial  Kuroshio 


with  sufficient  eastward  penetration  and  advection  of  heat 
as  well  as  cold  southward  currents  along  the  east  coast  of 
Japan  rather  than  the  warm  northward  flow  found  in  coarse- 
resolution  models.  Eddy-resolving  models  are  also  needed 
to  obtain  the  sharp  ocean  fronts  that  span  the  entire  North 
Pacific  [Hurlhurt  et  ai,  1996;  Hurlhurt  and  Metzger,  1998], 
In  the  North  Atlantic,  eddy-resolving  models  are  required  to 
accurately  simulate  the  Gulf  Stream  pathway  between  Cape 
Hatteras  and  the  Grand  Banks  and  to  obtain  the  associated 
large  nonlinear  recirculation  gyres  in  this  region  [Hurlhurt 
and  Hogan ,  2000;  Bryan  et  ai,  2007].  These  large  nonlinear 
gyres  change  the  large-scale  shape  of  the  subtropical  gyre 
into  a  C-shape.  East  of  the  Grand  Banks,  the  Gulf  Stream 
turns  northward  as  the  North  Atlantic  Current,  and  an  eddy- 
resolving  model  is  essential  for  the  model  to  advect  suffi¬ 
cient  heat  northward  and  then  eastward  across  the  Atlantic 
at  around  50°N  [Smith  et  ai,  2000].  These  are  examples  of 
regions  where  errors  caused  by  coarse  resolution  in  ocean 
global  circulation  models  (OGCMs)  would  lead  to  large  er¬ 
rors  in  climate  anomaly  or  climate  change  predictions  for 
North  America,  Europe,  the  Arctic,  and  Japan. 

1.2 .  Ocean  Weather  Prediction,  a  Brief  History  of  the 
First  10  Years 

GODAE  [Smith,  2000, 2006]  has  been  a  major  driving  force 
behind  a  multinational  effort  working  toward  the  develop¬ 
ment  of  eddy-resolving  global  ocean  prediction  capabilities 
in  different  countries.  Australia,  Britain,  France,  Japan,  and 
the  United  States  are  all  sponsoring  efforts  that  are  presently 
in  various  stages  toward  reaching  this  goal  (see  the  book  by 
Chassignet  and  Verron  [2006]  for  an  overview). 

The  Forecasting  Ocean  Assimilation  Model  (FOAM) 
system,  developed  by  the  Met  Office  UK,  became  the  first 
operational  global  ocean  prediction  system  in  1997,  but  at 
1°  resolution  it  was  non-eddy-resolving  [Bell  et  ai ,  2000]. 
Ten  years  later,  the  availability  of  sufficient  computer  re¬ 
sources  is  still  a  significant  limiting  factor  in  ocean  weather 
prediction.  Hence,  a  variety  of  approaches  and  trade-offs 
have  been  used  to  work  toward  the  eddy-resolving  global 
goal.  All  of  the  remaining  global-  and  basin-scale  predic¬ 
tion  systems  considered  here  are  either  eddy-permitting  or 
eddy-resolving,  and  all  include  some  form  of  altimeter  data 
assimilation.  SSH  from  satellite  altimetry'  is  the  key  avail¬ 
able  data  type  for  mapping  the  ocean  weather  in  deep  water 
[Hurlhurt,  1984]. 

Basin-scale  prediction  systems  have  been  a  common  step 
toward  the  global  goal,  either  nested  in  a  coarser  resolution 
global  system  (the  Met  Office  UK)  or  standalone  (the  re¬ 
mainder).  Such  systems  started  at  1/3°  resolution  in  Janu¬ 
ary  2001  (British  FOAM  Atlantic/Arctic  and  Indian  Ocean, 
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French  MERCATOR  North  and  tropical  Atlantic).  These 
were  soon  followed  by  higher  resolution  basin-scale  sys¬ 
tems  (1/9°  British  FOAM  North  Atlantic/Mediterranean 
Sea  since  2002,  operational  in  2004  [Bell  et  al .,  2006],  the 
French  MERCATOR  1/15°  North  Atlantic  plus  1/16°  Medi¬ 
terranean  Sea  operational  in  January  2003  [Bahurel,  2006], 
and  a  1/12°  Atlantic  system  run  in  demonstration  mode  by 
the  HYbrid  Coordinate  Ocean  Model  (HYCOM)  Consor¬ 
tium  in  the  United  States  since  July  2002  [C has  signet  et  al ., 
2007]). 

Some  of  the  basin-scale  systems  have  used  variable 
resolution  focused  on  a  region  of  interest.  The  Norwegian 
DIADEM  (later  TOPAZ)  Atlantic  and  Arctic  system  f Ber- 
lino  and  Evensen ,  2002]  has  been  running  since  2001  using 
a  curvilinear  grid  with  18-  to  35-km  resolution  (increased 
to  1 1-16  km  in  July  2007),  with  the  highest  resolution  fo¬ 
cused  on  the  northeast  Atlantic  and  Nordic  Seas.  Since  2001, 
the  Japan  Meteorological  Agency  has  been  running  a  daily 
operational  North  Pacific  assimilation-prediction  system 
COMPASS-K  (12°N  to  55°N)  with  variable  latitude  and 
longitude  resolution  as  high  as  1/4°  in  the  western  North  Pa¬ 
cific  around  Japan  [Kamachi  et  al .,  2004a,  2004b],  and  in 
2007,  the  Japanese  Meteorological  Research  Institute  began 
running  a  real-time  ocean  data  assimilation  and  prediction 
system  MOVE/M Rl.COM-NP  covering  the  Pacific  Ocean 
(15°S-65°N)  with  variable  latitude  and  longitude  resolution 
as  high  as  1/10°  near  Japan  [Usui  et  al .,  2006a,  2006b,  2008; 
Tsujino  et  al .,  2006].  The  US  National  Oceanic  and  Atmos¬ 
pheric  Administration/National  Centers  for  Environmental 
Prediction  (NOAA/NCEP)  has  been  running  a  real-time 
HYCOM  equatorial  and  North  Atlantic  prediction  system 
with  assimilation  of  satellite  altimeter  data  added  on  June  28, 
2007.  The  horizontal  resolution  is  about  4  to  18  km,  with  the 
highest  resolution  focused  in  the  Gulf  of  Mexico  and  along 
the  east  coast  of  North  America. 

Crosnier  and  Le  Provost  [2007]  performed  an  intercom¬ 
parison  study  of  five  basin-scale  prediction  systems  run  in 
five  different  countries  covering  the  North  Atlantic  and/or 
the  Mediterranean  Sea  [the  1/9°  British  FOAM  North  At- 
lantic/Mediterranean  Sea  system,  the  French  MERCATOR 
1/15°  North  Atlantic/1/16°  Mediterranean  Sea  system,  the  US 
1/12°  Atlantic  HYCOM  system,  the  HYCOM-based  Norwe¬ 
gian  TOPAZ  Atlantic/Arctic  system,  and  the  Italian  MFS  1/8° 
(later  1/16°)  Mediterranean  Sea  system].  Of  these,  three  use 
fixed-depth  r-level  models  and  two  use  HYCOM.  The  latter 
has  a  generalized  vertical  coordinate  that  allows  Fagrangian 
isopycnal  layers,  pressure  levels  (~r  levels),  terrain-following, 
and  other  vertical  coordinates  [ Bleck ,  2002;  Chassignet  et  al., 
2003;  Halliwell ,  2004;  Chassignet  et  al .,  2007]. 

Another  approach  is  that  of  the  Shallow- Water  Analysis 
and  Forecast  System  (SWAFS)  North  World  developed  at 


the  Naval  Oceanographic  Office  (NAVOCEANO).  SWAFS 
used  the  Princeton  Ocean  Model  with  terrain-following  co¬ 
ordinates  [Blutnberg  and  Mellor,  1987].  This  system  cov¬ 
ered  the  world  ocean  over  the  latitude  range  20°S  to  80°N 
with  -1/5°  (~23  km)  mid-latitude  resolution.  Effectively,  it 
was  a  convenient  way  to  run  multiple  basin-scale  models 
within  one  model  domain,  as  the  basins  are  disconnected 
at  the  southern  boundary  and  although  the  Indonesian  Ar¬ 
chipelago  is  included,  the  Indonesian  throughflow  is  largely 
blocked  because  of  the  southern  boundary.  The  SWAFS 
North  World  system  was  run  at  NAVOCEANO  from  1998 
to  2006,  but  never  became  operational.  A  SWAFS  system 
of  similar  design  for  the  Mediterranean  Sea  [Horton  et  al .. 
1997]  did  become  an  operational  prediction  system. 

Global  ocean  prediction  systems  that  make  various  com¬ 
promises  to  cope  with  limitations  in  available  comput¬ 
ing  power  have  also  been  developed.  A  1/4°  global  French 
MERCATOR  system  with  46  levels  in  the  vertical  has  been 
running  operationally  since  October  2005  [ Bahurel ,  2006], 
The  Australian  BFUEdink  global  ocean  prediction  system  has 
1/10°  resolution  on  a  90°  latitude  and  longitude  sector  sur¬ 
rounding  Australia,  variable  resolution  as  coarse  as  2°  else¬ 
where  and  47  levels  in  the  vertical  [Schiller  and  Smith,  2006]. 
The  BLUElink  system  became  operational  in  July  2007.  All 
of  the  preceding  systems  have  20  to  47  coordinate  surfaces  in 
the  vertical  direction. 

The  US  Navy  has  used  a  still  different  approach  in  its 
present  operational  global  ocean  prediction  systems,  a  linked 
two-model  approach  [Rhodes  et  al 2002].  One  model  has 
high  horizontal  resolution  but  low  vertical  resolution,  while 
the  other  has  high  vertical  resolution  but  lower  horizontal 
resolution  The  earliest  version  of  this  system  used  only 
the  model  with  low  vertical  resolution,  the  Naval  Research 
Faboratory  (NRL)  Layered  Ocean  Model  (NLOM).  It  had 
six  Lagrangian  layers  in  the  vertical  and  1/4°  resolution  glo¬ 
bally,  excluding  the  Arctic  and  most  shallow  water  [Metzger 
et  al .,  1998a,  1998b].  It  was  the  first  global  ocean  weather 
prediction  model  to  assimilate  real-time  satellite  altimeter 
data  (from  TOPEX/Poseidon  and  ERS-2)  and  ran  daily  in 
real  time  from  1997  to  2000  at  the  Fleet  Numerical  Mete¬ 
orology  and  Oceanography  Center  (FNMOC)  but  never  be¬ 
came  operational  In  2000,  NAVOCEANO  became  the  US 
Navy  center  for  numerical  prediction  of  the  ocean  weather. 
A  mixed  layer  (seventh  layer)  [Waller aft  et  al 2003]  was 
added  to  the  NL  OM-based  component  of  a  system  for 
NAVOCEANO,  which  began  running  daily  in  real  time  with 
1/16°  resolution  on  October  18,  2000.  It  ran  operationally 
from  September  27,  2001  to  March  11,  2006,  thereby  be¬ 
coming  the  first  operational  eddy- resolving  global  ocean 
prediction  system  [Smedstad et  al.,  2003].  A  1/32°  version  of 
this  system  began  running  in  near  real  time  on  November  1, 
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2003,  real  time  since  March  1, 2005,  and  replaced  the  1/16° 
system  as  the  operational  system  on  March  6,  2006  [Sh river 
et  ai ,  2007]. 

The  second  component  of  the  NAVOCEANO  system 
uses  the  Navy  Coastal  Ocean  Model  (NCOM)  and  has  1/8° 
(~15  km)  mid-latitude  resolution,  1/6°  (~19  km)  tropical  res¬ 
olution,  40  levels  in  the  vertical,  and  is  fully  global  [ Barron 
et  al. ,  2006;  Kara  et  ai ,  2006].  It  has  been  running  daily  in 
real  time  since  October  2001  and  became  operational  on  Feb¬ 
ruary  19,  2006.  The  NLOM  component  is  used  to  assimilate 
SSH  data  along  altimeter  tracks  with  a  model  forecast  as 
the  first  guess  for  each  analysis  cycle  and  it  is  used  to  make 
30-day  ocean  weather  forecasts.  It  also  assimilates  model- 
independent  analyses  of  sea  surface  temperature  (SST).  The 
NCOM  component  assimilates  steric  SSH  anomalies  from 
NLOM  in  the  form  of  synthetic  temperature  and  salinity  (T  & 
S)  profiles  [Rhodes  et  al .,  2002].  These  profiles  are  obtained 
from  surface  dynamic  height  and  SST  anomalies  by  means 
of  regression  equations  for  subsurface  temperature  and  T  & 
S  relations  for  salinity,  both  using  statistics  derived  from  the 
historical  hydrographic  data  base  [ Carnes  et  al .,  1990;  Fox  et 
al. ,  2002].  This  two-model  global  system  requires  less  com¬ 
puter  power  than  a  single  global  model  with  both  high  hori¬ 
zontal  and  high  vertical  resolution,  which  is  the  goal  of  most 
GODAE  participants. 

The  US  Navy  began  running  the  first  eddy-resolving  glo¬ 
bal  system  with  high  vertical  resolution  in  near  real  time  on 
December  22,  2006,  and  it  has  been  running  in  real  time 
since  February  16,  2007.  This  system  uses  a  fully  global 
configuration  of  HYCOM  with  1/12°  equatorial  (~7  km  mid¬ 
latitude)  resolution  and  32  layers  in  the  vertical.  Each  day, 
it  performs  a  5-day  hindcast  (assimilation  of  data  before  the 
real-time  data  window)  at  1-day  increments  to  pick  up  de¬ 
layed  data,  and  it  makes  a  five-day  forecast.  In  addition,  tests 
of  30-day  forecasts  have  been  performed. 

Although  limitations  on  ocean  model  domain  size  and 
resolution  have  been  essential  in  dealing  with  limitations 
in  available  computer  power,  computing  requirements  have 
also  been  economized  through  the  choice  of  data  assimila¬ 
tion  techniques  and  assimilated  data  sets.  In  addition,  while 
most  of  the  preceding  systems  run  daily  in  real  time,  com¬ 
puting  requirements  for  other  systems  have  been  reduced  by 
running  in  near  real  time  and  updating  once  a  week,  rather 
than  daily. 

1.3.  Ocean  Prediction  in  Relation  to  Classes  of  Ocean 
Response  to  Atmospheric  Forcing 

Table  1  is  a  summary  of  some  useful  classes  of  ocean  re¬ 
sponse  to  atmospheric  forcing  that  cover  many  phenomena  of 
potential  interest  in  ocean  prediction  In  developing  an  ocean 


prediction  system,  it  is  essential  to  consider  these  classes  and 
their  implications  for  ( 1 )  data  requirements,  (2)  the  timescales 
of  oceanic  predictive  skill,  and  (3)  prediction  system  design. 

In  the  case  of  class  1  phenomena,  characterized  by  strong 
and  rapid  development  (less  than  a  week),  ocean  model 
simulation  skill  plays  a  critical  role  in  converting  atmos¬ 
pheric  forcing  into  oceanic  information  The  timescale  for 
oceanic  forecast  skill  is  largely  determined  by  the  timescale 
for  atmospheric  predictive  skill,  and  the  influence  of  anom¬ 
alies  in  the  initial  state  rapidly  decreases.  The  timescale  for 
oceanic  predictive  skill  may  exceed  that  for  the  atmosphere 
to  some  extent  because  of  lagged  responses,  the  influence 
of  geometric  or  topographic  constraints,  or  the  initiation  of 
class  4  free  waves  or  ocean  eddies.  Such  phenomena  are 
commonly  found  in  shallow  water  where  rapid  barotropic 
responses  are  particularly  important,  within  the  surface 
mixed  layer  or  Ekman  layer,  in  equatorial  and  coastal  wave 
guides,  and  in  the  vicinity  of  hurricanes  or  other  strong  wind 
events  (e.g.,  wind-driven  eddies  in  the  Gulf  of  Tehuante¬ 
pec).  In  section  2.7,  some  examples  are  given  in  which 
class  1  anomalies  lead  to  class  4  anomalies  or  where  class 
1  and  class  4  phenomena  interact  to  generate  new  class  4 
features. 

Nowcasting  and  forecasting  of  class  1  phenomena  place 
high  demands  on  the  accuracy  of  the  ocean  model  and  the  at¬ 
mospheric  forcing.  In  shallow  water  and  near  coastal  bounda¬ 
ries,  there  are  additional  demands  for  accurate  topography 
and  coastlines,  tides,  and  fine-resolution  atmospheric  forc¬ 
ing.  Near  coastlines,  it  is  important  that  atmospheric  values 
over  land  are  not  used  as  forcing  at  ocean  model  grid  points, 
a  common  problem  in  fine-resolution  ocean  models  forced  by 
a  coarser  resolution  atmospheric  product  [Kara  et  al. ,  2007]. 
When  ocean  data  are  assimilated,  the  model  time  and  the  ob¬ 
servation  time  must  closely  match  because  of  the  rapid  evolu¬ 
tion;  but  even  so,  the  ocean  data  may  rapidly  lose  impact  after 
assimilation. 

Class  2  phenomena  (slower  and  indirectly  forced)  tend  to 
be  a  largely  nondeterm inistic  response  to  atmospheric  forc¬ 
ing  due  to  flow  instabilities,  for  example,  the  generation  of 
mesoscale  eddies  and  the  meandering  of  ocean  currents  and 
fronts.  Thus,  ocean  data  assimilation  is  essential  to  nowcast 
and  forecast  these  features,  and  the  nowcast  and  forecast 
results  are  relatively  insensitive  to  the  atmospheric  forcing, 
allowing  forecasts  up  to  a  month  or  more.  SSH  anomalies 
measured  by  satellite  altimetry  are  the  key  data  that  all  of  the 
GODAE  participants  use  to  map  these  mesoscale  features, 
as  well  as  many  other  features,  including  El  Nino  and  La 
Nina  events.  SSH  is  a  particularly  useful  data  type  observable 
from  space  because  (1 )  it  is  geostrophically  related  to  surface 
currents,  (2)  it  is  an  integral  measure  of  subsurface  variations 
in  T  &  S  (e.g.,  the  depth  of  the  thermocline)  plus  a  bottom 
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Table  1.  Classes  of  Ocean  Response  lo  Atmospheric  Forcing.3 

Class  Examples  Implications 


1.  Strong,  rapid  (<l  week), 
and  direct 


2.  Slower  (weeks  to  months) 
and  indirect 


Surface  waves,  storm  surges,  and  rapid 
variations  within  the  upper  mixed  layer, 
beneath  hurricanes,  in  coastal  and  equatorial 
upwelling,  and  often  the  onset  of  equatorial 
and  coastal  trapped  waves 
Mesoseale  eddies,  meandering  currents, 
frontal  locations,  features  related  to  flow 
instabilities  on  the  mesoseale 


3.  Slow  (weeks  to  years)  El  Nino;  much  of  the  tropical  ocean 

direct  integrated  response  circulation,  gyres,  persistent  features  associated 

with  geometric  and  topographic  constraints 


4.  Free  propagation  of  existing  Equatorial  and  coastal  I  y  trapped  waves,  Rossby 
features  (weeks  to  years)  waves  and  stable  isolated  eddies  generated 

under  Classes  1 , 2,  and  3 


Forecasts  are  short  range;  limited  by 
atmospheric’  predictive  skill  Less  sensitive 
to  errors  in  the  initial  state:  more  sensitive  to 
errors  in  forcing 

Forecast  skill  up  to  a  month  or  more,  more 
sensitive  to  errors  in  the  initial  state:  less 
sensitive  to  errors  in  forcing,  ocean  data 
assimilation  is  essential;  variability  statistics 
may  be  predicted  via  simulation 
Long  range  forecasts  possible;  sensitive  only 
to  errors  in  forcing  on  long  time  scales;  now- 
casting  and  forecasting  feasible  using  ocean 
models  with  sparse  ocean  data 
Reduced  sensitivity  to  onset  mechanisms  and 
extended  forecast  skill,  weeks  to  years 


“This  table  summarizes  the  classes  of  ocean  response  lo  atmospheric  forcing  most  relevant  to  ocean  climate  and/or  ocean  weather  predie 
lion,  the  first  three  in  order  of  increasing  time  scale.  The  fourth  class  is  initiated  by  anomalies  generated  under  the  preceding  three  classes. 
Class  2  is  a  close  analog  of  atmospheric  cyclones,  anticyclones,  and  jet  streams,  but  with  larger  timescales  and  much  smaller  space  scales, 
the  latter  largely  due  to  the  relative  atmospheric  and  oceanic  first  internal  radii  of  deformation.  A  fifth  class  includes  tides  and  tsunamis 
that  are  (mostly)  not  in  response  to  the  atmosphere.  Other  classes  are  not  considered  here,  for  example,  most  inertia-gravity  waves  and 
line-scale  flow  instabilities  that  may  result  from  atmospheric’  forcing.  Adapted  from  Hurlburt  [1984], 


pressure  anomaly,  and  (3)  the  oceanic  vertical  structure  is 
generally  low  mode  Thus,  it  has  been  possible  to  develop  a 
number  of  techniques  to  project  the  SSH  anomalies  down¬ 
ward,  in  some  cases  including  the  bottom  pressure  anomaly 
(the  nonsteric  contribution  to  SSH).  These  techniques  range 
from  model  dynamics  [Hurlburt,  1986],  to  model  statistics 
[Hurlburt  et  al. ,  1990],  lo  a  potential  vorticily  constraint 
[Cooper  and  Haines,  1996],  to  data  assimilation  covariances, 
to  synthetic  T  &  S  profiles  estimated  from  a  combination  of 
SSH  and  SST  related  by  regression  to  the  historical  hydro- 
graphic  data  base  [DeMey  and  Robinson ,  1 987;  Carnes  et  al ., 
1 990;  Fox  et  al.,  2002]. 

The  regions  of  the  global  ocean  dominated  by  class  2  van- 
ability  can  be  estimated  by  calculating  the  fraction  of  variabil¬ 
ity  of  different  ocean  variables  [c.g.,  SSH,  SST,  eddy  kinetic 
energy  (EKE),  and  bottom  pressure]  that  is  a  nondeterministic 
response  to  atmospheric  forcing  in  eddy-resolving  model  sim¬ 
ulations  [Metzger  et  al .,  1994;  Hurlburt  and  Metzger,  1998; 
Metzger  and  Hurlburt ,  200 1 ;  Melsom  et  a l .,  2003;  Hogan  and 
Hurlburt ,  2005].  Although  areas  that  are  highly  nondetermin¬ 
istic  are  not  necessarily  regions  of  high  SSH  variability  or 
EKE,  it  is  relevant  to  evaluate  the  model  against  its  ability  to 
simulate  regions  of  high  variability.  In  addition,  ocean  fore¬ 
cast  skill  beyond  the  range  of  atmospheric  predictive  skill 
can  be  verified.  The  sensitivity  of  ocean  forecast  skill  to  at¬ 


mospheric  forcing  can  be  tested  versus  persistence  (a  fore¬ 
cast  of  no  change)  and  by  comparing  the  skill  of  forecasts 
performed  using  analysis  quality  forcing  versus  forecasts  re¬ 
verting  toward  climatological  forcing  beyond  the  end  of  the 
forecast  forcing  [Smedstad  et  al .,  2003;  Shriver  et  al .,  2007, 
section  2  of  this  chapter]  (relevant  to  all  four  classes  of  oce¬ 
anic  response). 

Class  3  phenomena  (slow,  direct  integrated  responses  to 
atmospheric  forcing  on  timescales  of  weeks  to  years)  arc 
perhaps  more  relevant  to  seasonal  to  interannual  prediction 
(e.g.,  El  Nino  and  La  Nina  events)  using  coupled  ocean 
atmosphere  models.  However,  class  3  phenomena  are  also 
relevant  to  uncoupled  models  and  shorter  range  predictions. 
They  are  relevant  to  both  types  of  prediction  when  consider¬ 
ing  the  role  of  uncoupled  models  in  determining  the  mean 
SSH  added  to  SSH  anomalies  from  satellite  altimetry,  a  topic 
discussed  in  section  2.6  The  required  mean  SSH  contains  a 
wide  range  of  space  scales,  some  of  which  are  adequately  re¬ 
solved  only  in  simulations  by  a  high-resolution  model.  This 
places  a  heavy  burden  on  the  accuracy  of  the  ocean  model 
and  the  atmospheric  forcing. 

In  addition,  class  3  phenomena  are  important  in  nowcast¬ 
ing  and  forecasting  some  ocean  features  that  are  inadequately 
observed,  for  example,  for  an  eddy-resolving  model,  the 
seasonal  cycle  of  mixed  layer  depth  (MI  D)  and  anomalies 
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in  MLD  lasting  longer  than  the  class  1  timescale  of  skill¬ 
ful  atmospheric  weather  prediction,  examples  where  a  class 
3  response  to  atmospheric  forcing  will  have  a  much  larger 
impact  than  assimilation  of  ocean  data.  Presently,  -3000 
Argo  floats  [Roemmich  et  al .,  2004]  profile  T  &  S  down  to 
typical  depths  of  2000  m  with  high  vertical  resolution.  These 
floats  constitute  an  accurate  observing  system  for  measuring 
subsurface  T  &  S,  including  MLD,  but  one  that,  while  suit¬ 
able  for  climate,  is  too  sparse  horizontally  and  temporally  to 
constrain  an  eddy-resolving  ocean  model  and  too  sparse  to 
constrain  MLDs  in  any  model  on  class  I  timescales.  As  dis¬ 
cussed  by  Hurlburt  et  al.  [2001],  3000  Argo  floats,  arrayed 
like  altimeter  tracks,  could  provide  horizontal  coverage  at  a 
single  level  (e.g.,  MLD  or  surface  dynamic  height)  which 
is  only  -4%  of  that  from  one  Jason- 1  altimeter  over  its  10- 
day  orbital  repeat  cycle.  This  comparison  is  based  on  plac¬ 
ing  Argo  floats  along  over-water  altimeter  tracks  at  50-km 
intervals  and  profiling  at  10-day  intervals. 

Class  4  phenomena  consist  of  free  propagation  of  ocean 
features  that  exist  in  the  initial  state  of  the  forecast,  that  is, 
features  generated  earlier  as  class  1,  2,  or  3  phenomena  In 
addition  to  examples  mentioned  in  the  discussion  of  class  1 
phenomena,  Kelvin  waves  generated  by  El  Ninos  have  prop¬ 
agated  cyclonically  around  the  North  Pacific  as  far  as  Alaska 
and  the  Kamchatka  Peninsula  of  Russia  in  both  observations 
and  models  [Metzger  et  al .,  1998a;  Melsom  et  al .,  2003]. 
Jacobs  et  al.  [1994]  even  demonstrated  a  decade-long  fore¬ 
cast  of  a  trans-Pacific  nondispersive  internal  Rossby  wave, 
generated  from  the  1982-1983  El  Nino,  which  eventually 
displaced  the  Kuroshio  Extension  east  of  Japan  The  fore¬ 
cast  used  climatological  atmospheric  forcing  (a  correspond¬ 
ing  simulation  with  interannual  forcing  was  also  run),  and 
the  forecast  was  verified  using  satellite  altimetry  and  other 
data.  In  addition  to  wind-generated  eddies  (like  Tehuantepec 
eddies),  stable  isolated  eddies  that  are  products  of  flow  insta¬ 
bility,  but  which  escape  any  region  of  strong  flow  instability, 
may  belong  to  class  4.  Such  eddies  can  last  for  over  a  year 
[Lai  and  Richardson ,  1 977],  although  they  are  unlikely  to  be 
predictable  for  nearly  that  long. 

Nowcasting  is  an  essential  step  in  ocean  monitoring  and 
prediction  of  all  the  classes  of  oceanic  response  in  Table  1 . 
To  perform  a  nowcast,  a  numerical  model  forecast  is  used 
as  the  first  guess  for  assimilation  of  new  data.  This  use  of 
the  forecast  allows  the  model  ( 1 )  to  help  fill  in  the  temporal 
and  spatial  gaps  between  observations  (including  tempo¬ 
ral  gaps  between  the  repeat  cycles  of  satellite  altimetry) 
and  make  effective  use  of  delayed  data  by  exploiting  its 
predictive  skill,  (2)  to  help  convert  the  better  observed  sur¬ 
face  fields  into  subsurface  structure  [Hurlburt,  1986],  (3) 
to  convert  the  better  observed  atmospheric  forcing  func¬ 
tions  into  useful  oceanic  responses,  and  (4)  to  apply  bottom 


topography,  coastline  geometry,  and  the  ocean  surface  as 
physical  and  dynamical  constraints  with  associated  bound¬ 
ary  layers.  Partially  in  conjunction  with  a  barotropic  model 
[Carrere  and  Lyard ,  2003],  the  model  component  of  the 
ocean  prediction  system  can  also  help  separate  the  steric 
and  nonsteric  contributions  to  SSH,  a  significant  capabil¬ 
ity  in  the  assimilation  of  satellite  altimeter  data  [Rhodes  et 
al .,  2002].  Examples  of  the  preceding  model  contributions 
to  nowcast  and  forecast  skill  (without  a  reference)  are  il¬ 
lustrated  in  section  2. 

There  is  an  important  difference  between  data  assimila¬ 
tion  in  oceanography  and  meteorology.  In  oceanography, 
there  is  a  greater  burden  on  the  models  to  extract  informa¬ 
tion  from  the  data,  including  use  of  the  models  to  extract 
oceanic  responses  to  atmospheric  forcing  as  well  as  con¬ 
verting  surface  information  into  subsurface  information, 
which  is  sparsely  observed.  In  addition,  an  eddy-resolving 
global  ocean  model  should  play  a  major  role  in  determin¬ 
ing  the  mean  SSH  used  in  assimilation  of  satellite  altimeter 
data  (discussed  in  section  2),  and  it  may  be  used  as  a  source 
of  covariances  and  other  statistics  for  data  assimilation. 
Therefore,  models  play  a  crucial  role  in  making  effective 
use  of  the  data  sources  and  in  reducing  data  requirements  in 
ocean  nowcasting  in  addition  to  their  essential  role  in  ocean 
forecasting. 

Other  applications  of  eddy-resolving  ocean  models  in  ocean 
prediction  include  understanding  model  and  ocean  dynam¬ 
ics,  assessing  and  improving  the  realism  of  the  ocean  model, 
studying  ocean  model  resolution  requirements,  performing 
oceanic  predictability  studies,  investigating  the  classes  of  re¬ 
sponse  to  atmospheric  forcing  in  different  dynamical  regimes, 
and  providing  realistic  simulated  data  (a)  for  data  assimila¬ 
tion  studies,  (b)  for  studies  of  data  base  requirements  and  ob¬ 
serving  system  simulation  experiments,  and  (c)  for  spinning 
up  models  to  the  initial  time  when  data  assimilation  begins. 
Eddy-resolving  models  can  also  be  used  to  investigate  model 
sensitivity  to  the  choice  of  atmospheric  forcing  product,  al¬ 
though  a  linear  barotropic  numerical  model  with  realistic  deep 
water  boundaries  (essentially  Sverdrup  [1947]  interior  flow 
with  Munk  [1950]  viscous  western  boundary  layers  consist¬ 
ent  with  the  Godfrey  [1989]  island  rule)  can  be  very  useful 
for  this  purpose  with  much  lower  computational  require¬ 
ments  [Townsend  et  al .,  2000;  Hogan  and  Hurlburt ,  2005]. 
A  number  of  the  issues  raised  in  this  subsection  are  illustrated 
and  discussed  in  section  2. 

2.  EDDY-RESOLVING  GLOBAL  OCEAN 
PREDICTION  SYSTEMS 

Here,  we  focus  on  the  only  two  existing  eddy-resolving 
global  ocean  prediction  systems  to  date,  the  1/16°  (later 
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1/32°)  NLOM-based  system  with  low  vertical  resolution 
and  the  more  reeent  1/12°  global  HYCOM-based  system, 
which  has  both  high  horizontal  and  high  vertical  resolution. 
HYCOM  is  fully  global,  unlike  NLOM,  whieh  excludes  the 
Aretie  and  most  shallow  water.  The  results  presented  are 
used  in  the  illustration  and  diseussion  of  some  key  issues  in 
oeean  weather  prediction. 

2. 1.  The  HYCOM  and  NLOM  Prediction  Systems 

Both  HYCOM  and  NLOM  have  a  free  surfaee  and  use  a 
Lagrangian  vertical  coordinate  to  model  the  oeean  interior. 
However,  simplifications  are  made  in  NLOM  that  make  it 
much  more  computationally  efficient,  but  which  limit  its 
range  of  application  and  versatility,  while  the  generalized 
vertical  coordinate  and  the  allowance  of  zero  thiekness  lay¬ 
ers  make  HYCOM  extremely  versatile  and  suitable  for  a 
w  ide  range  of  oeean  modeling  applications. 

General  details  of  NLOM  design  and  progressive  devel¬ 
opment  are  discussed  by  Hurlhurt  and  Thompson  [1980], 
Wall  craft  [1991],  Wall  craft  and  Moore  [1997],  Moore  and 
W allcraft  [  1 998],  and  Waller  aft  et  al.  [2003].  The  last  of  these 
discusses  the  development  of  NLOM  as  a  thermodynamic 
model  with  a  Kraus  and  Turner  [1967]-type  bulk  mixed 
layer  and  SST.  This  version  of  NLOM  was  used  in  the  1/16° 
and  1/32°  global  prediction  systems  listed  in  section  1.2. 
More  specifically,  the  resolution  for  eaeh  variable  is  1/16° 
(1/32°)  in  latitude  and  45/512°  (45/1024°)  in  longitude  or 
~7  km  (3^4  km)  at  mid-latitudes  in  the  1/16°  (1/32°)  system. 
The  NLOM  domain  is  nearly  global,  extending  from  72°S 
to  65°N,  with  the  model  boundary  generally  following  the 
200-m  isobath  with  a  few  exceptions,  like  the  shallow  straits 
connecting  the  Sea  of  Japan  (East  Sea)  to  the  Pacific  Ocean 
NLOM  uses  vertically  eompressed  but  otherwise  realistic 
bottom  topography  confined  to  the  lowest  layer  [Hurlburt  et 
al..  Pathways  of  upper  oeean  currents  and  fronts:  Steering  by 
the  topographical ly-eonstrained  abyssal  circulation  and  the 
role  of  flow  instabilities,  applicability  of  a  two-layer  theory 
to  low  versus  high  vertical  resolution  models,  submitted  to 
Dynamics  of  the  Atmospheres  and  Oceans ,  2008  hereinafter 
referred  to  as  Hurlburt  et  al.,  submitted  manuscript,  2008]. 
Sill  depths  of  straits  that  are  shallower  than  the  bottom  layer 
are  maintained  by  constraining  the  flow  to  small  values  be¬ 
low  the  sill  depth  [Metzger  and  Hurlburt ,  1 996].  The  redueed 
amplitude  of  the  topography  generally  does  not  have  a  sig¬ 
nificant  impact  on  the  abyssal  current  steering  of  upper  ocean 
current  pathways  (an  essential  capability  of  eddy-resolving 
oeean  models  which  simulate  eddy-driven  abyssal  currents) 
because  even  low-amplitude  topographic  features  ean  eon- 
strain  the  pathways  of  abyssal  flow  [ Hurlburt  and  Metzger , 
1998;  Hurlburt  et  al.,  submitted  manuscript,  2008].  Deep  wa¬ 


ter  formation  outside  the  model  domain  is  included  via  pre¬ 
scribed  observational ly  based  transports  through  four  straits 
in  the  northern  boundary  of  the  North  Atlantic,  northw  ard  in 
the  upper  layers  and  southward  in  the  bottom  layer  [Shriver 
and  Hurlburt,  1997]. 

The  history'  of  NLOM  use  in  ocean  prediction  systems  is 
diseussed  in  seetion  1 .2.  Eaeh  day,  the  NLOM  prediction  sys¬ 
tem  performs  a  hindeast  at  daily  increments  to  piek  up  de¬ 
layed  altimeter  data  (a  5-day  hindeast  for  the  1/32°  system,  a 

3- day  hindeast  for  the  former  1/16°  system),  and  it  makes  a 

4- day  forecast  (plus  a  30-day  foreeast  onee  a  week)  [Smed- 
stad  et  al .,  2003;  Shriver  et  al .,  2007].  NLOM  assimilates 
along-track  satellite  altimeter  data  using  the  model  as  a  first 
guess  for  the  SSH  analysis,  and  it  assimilates  SST  in  the 
form  of  daily  operational  model-independent  SST  analyses 
from  N  AVOCEANO  [Barron  and  Kara ,  2006].  The  SSH  as¬ 
similation  proeess  consists  of  an  optimum  interpolation  (Ol) 
deviation  analysis  from  the  model  first  guess  and  an  empiri¬ 
cal  orthogonal  function  regression  technique  based  on  model 
statistics  to  project  the  SSH  updates  downward,  including 
to  the  abyssal  layer  [Hurlburt  et  al .,  1990].  These  updates 
arc  gcostrophically  balanced  outside  the  equatorial  region 
and  inserted  incrementally  to  further  reduce  inertia-gravity 
wave  generation.  Anisotropie,  spatially  vary  ing  mesoseale 
eovariance  functions  detenu ined  from  TOPEX/Poseidon 
and  ERS-2  altimeter  data  [.Jacobs  et  al .,  2001]  are  used  in 
the  Ol  analysis.  The  SST  assimilation  consists  of  relaxing 
the  analysis  into  the  model. 

HYCOM  has  been  developed  as  a  collaborative  multi- 
institutional  effort  starting  from  the  Miami  Isopycnie  Coordi¬ 
nate  Ocean  Model  (MICOM)  using  the  theoretical  foundation 
set  forth  by  Blech  and  Boudra  [1981],  Blech  and  Benjamin 
[1993],  and  Blech  [2002].  HYCOM  uses  a  generalized  verti¬ 
cal  coordinate  that  is  normally  isopyenal  in  the  open  strati¬ 
fied  ocean,  with  a  dynamically  smooth  transition  to  pressure 
coordinates  (~z-levels)  in  the  unstratified  mixed  layer  and  to 
a  (terrain-following)  coordinates  in  shallow  water,  but  it  is 
not  limited  to  these  types.  The  transition  between  coordinate 
types  is  made  dynamically  in  spaee  and  time  using  the  lay  ered 
continuity  equation  and  a  hybrid  coordinate  generator.  This 
generalized  coordinate  approach  retains  particular  advan¬ 
tages  associated  with  the  different  vertical  coordinate  types: 
( 1 )  the  retention  of  water  mass  charaeteristies  for  eenturies  (a 
eharaeteristie  of  the  isopyenal  coordinate  interior),  (2)  high 
vertical  resolution  in  the  surface  mixed  layer  and  unstratified 
or  weakly  stratified  regions  of  the  ocean  (a  eharaeteristie  of 
r-level  coordinates),  and  (3)  high  vertical  resolution  in 
coastal  regions  (a  eharaeteristie  of  terrain-following  coordi¬ 
nates)  [Chassignet  et  al .,  2003].  The  generalized  coordinate 
also  facilitates  accurate  transition  between  deep  and  shallow 
water  and  allows  the  use  of  a  sophisticated  embedded  mixed- 
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layer  model  in  an  ocean  model  with  an  isopycnal  interior,  for 
example,  a  K-profile  parameterization  [Large  et  al .,  1997]  in 
the  ocean  prediction  system  For  the  1/12°  global  HYCOM, 
the  model  resolution  south  of  47°N  is  exactly  0.08°  by  0.08 
cos(latitude)0  for  each  variable  (~7  km  at  mid-latitudes)  and 
is  described  by  the  approximate  equatorial  resolution.  A  bi¬ 
polar  grid  is  used  north  of  47°N,  resulting  in  a  tripole  global 
grid  [ Murray ,  1996],  which  gives  ~3.5-km  resolution  in  the 
Arctic.  The  model  includes  32  hybrid  coordinate  surfaces  in 
the  vertical  direction  and  is  run  with  thermobancity  [Sun  et 
al..  1999;  Chassignet  et  al .,  2003;  Hallberg. ,  2005]  using  a 
reference  depth  of  2000  db  for  potential  density,  a  model  con¬ 
figuration  known  as  a2*. 

The  1/12°  global  HYCOM  prediction  system  has  been 
running  daily  in  near  real  time  since  December  22,  2006  and 
in  real  time  since  February  16,  2007.  Even  more  than  for 
the  NLOM  system,  because  HYCOM  includes  shallow  wa¬ 
ter,  doubling  the  resolution  is  desirable  when  the  computing 
power  becomes  available  The  impact  of  model  resolution 
is  a  recurring  theme  in  the  rest  of  section  2.  The  data  are  as¬ 
similated  into  HYCOM  using  the  Navy  Coupled  Ocean  Data 
Assimilation  (NCODA)  system  of  Cummings  [2005]  with 
a  modification  described  in  section  2.8.  Within  NCODA,  a 
HYCOM  forecast  is  used  as  the  first  guess  for  a  multivari¬ 
ate  optimum  interpolation  (MVOl)  [Daley,  1991]  analysis 
of  SSH  track  data  from  satellite  altimetry,  SST,  and  T  & 
S  profiles.  The  analysis  is  then  inserted  incrementally  into 
HYCOM.  In  the  analysis,  the  method  of  Cooper  and  Haines 
[1996]  is  used  for  downward  projection  of  SSH  below  the 
mixed  layer,  although  downward  projection  of  SSH  and  SST 
using  synthetic  T  &  S  profiles  [Fox  et  al.,  2002]  is  also  an 
option  in  NCODA.  The  updates  to  the  mass  field  are  geos- 
trophically  balanced  outside  an  equatorial  band. 

2.2.  1/16°  Eddy- Resolving  Versus  1/4°  Eddy- Permitting 
Models  in  Nowcasting  and  Forecasting  the  Kuroshio 
Extension 

Plate  1  illustrates  the  relative  dynamical  interpolation 
(nowcast)  and  forecast  skill  of  1/4°  eddy-permitting  ver¬ 
sus  1/16°  eddy-resolving  six-layer  versions  of  NLOM  in 
the  Kuroshio  Extension  region  east  of  Japan,  as  verified  by 
NAVOCEANO  operational  1/8°  Modular  Ocean  Data  As¬ 
similation  System  (MODAS)  model-independent  SST  analy¬ 
ses  [Barron  and  Kara ,  2006]  with  the  color  bar  designed  to 
highlight  the  Kuroshio  pathway.  The  1/16°  Pacific  prediction 
system  discussed  here  is  the  forerunner  of  the  1/16°  global 
NLOM  prediction  system  which  later  became  operational  at 
NAVOCEANO,  except  that  it  lacked  a  mixed  layer  and  SST 
[Hurlburt  et  al. ,  2000;  Metzger  et  al.,  2000].  The  1/4°  global 
system  is  one  of  two  versions  run  in  real  time  at  FNMOC 


(section  1.2).  Results  for  January  1  and  January  15,  1999  are 
shown  for  each  product.  Two- week  forecasts  from  January 
I  for  each  model  can  be  compared  with  the  three  analyses 
for  January  1 5.  The  model  nowcasts  were  updated  daily  with 
assimilation  of  altimeter  track  data  from  TOPEX/Poseidon 
and  ERS-2  using  a  3-day  window  (red  tracks  overlaid  on 
two  of  the  plots)  and  the  methodology  outlined  in  section  2. 1 
for  NLOM.  Because  the  geoid  of  the  Earth  is  not  adequately 
known,  only  the  altimetric  deviations  from  their  own  mean 
are  useful  in  the  assimilation.  In  the  model  assimilations,  a 
slightly  modified  1993  1997  mean  from  the  corresponding 
model  was  added  to  these  deviations.  Generally,  the  two 
means  agree  closely,  but  the  1/16°  model  gives  a  sharper  de¬ 
piction  of  mean  currents  and  a  sharper  depiction  of  the  Kuro¬ 
shio  in  the  nowcast  and  forecast  results. 

Several  features  are  clearly  evident  in  the  SST  and  the 
model  results.  Here,  we  focus  on  the  sharp  meander  between 
155°E  and  160°E  on  January  1,  1999,  which  evolved  into 
a  distinctive  eddy  detachment  in  progress  on  January  15. 
This  sequence  is  captured  with  striking  agreement  between 
the  SST  analyses,  the  1/16°  model  nowcasts,  and  the  1/16° 
model  forecast.  The  1/4°  model  represents  the  sharp  mean¬ 
der  on  January  1  quite  well,  but  the  1/4°  nowcast  for  the  de¬ 
taching  eddy  on  January  15  is  poor  and  the  January  15  1/4° 
model  forecast  for  this  feature  is  a  very  unrealistic  meander 
contraction. 

T  hese  results  are  a  demonstration  of  ocean  model  eddy¬ 
resolving  nowcast  and  forecast  skill  using  satellite  altimeter 
data.  In  particular,  the  results  demonstrate  (1)  that  satellite 
altimetry  can  be  an  effective  observing  system  for  map¬ 
ping  and  forecasting  mesoscale  ocean  features,  (2)  that  an 
ocean  model  with  high  enough  horizontal  resolution  can  be 
a  skillful  dynamical  interpolator  of  satellite  altimeter  data 
in  depicting  mesoscale  oceanic  variability,  and  (3)  that  an 
eddy-resolving  ocean  model  can  provide  skillful  forecasts 
of  mesoscale  variability  when  a  model  with  assimilation  of 
altimeter  data  is  used  to  define  the  initial  state,  in  this  exam¬ 
ple  for  the  duration  of  the  31 -day  forecast  [Hurlburt  et  al., 
2000].  It  is  worth  noting  that  this  was  the  very  first  forecast 
performed  with  the  1/16°  Pacific  model  used  here,  not  the 
best  of  many.  In  contrast,  the  1/4°  global  model  had  insuf¬ 
ficient  horizontal  resolution,  lower  skill  as  a  dynamical  in¬ 
terpolator,  and  poor  forecast  skill,  with  a  strong  tendency  to 
damp  the  variability.  The  broadening  of  the  Kuroshio  in  the 
1/4°  model  means  the  current  speeds  are  lower,  and  dynami¬ 
cally,  the  current  is  not  as  inertial,  a  clear  indication  that  one 
should  anticipate  a  substantial  impact  on  model  forecast  skill 
and  skill  as  a  dynamical  interpolator. 

The  1/16°  Pacific  model  used  for  assimilation  of  real  al¬ 
timeter  data  in  Plate  1  was  also  used  in  testing  assimilation  of 
simulated  altimeter  track  data  derived  from  the  model  along 
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ERS-2,  GFO,  and  both  sets  of  TOPEX/Poseidon  repeat 
tracks.  Simulated  altimeter  data  from  model  year  1994  were 
assimilated  into  the  model  starting  from  an  initial  state  in 
model  year  1997,  a  time  when  the  model  Kuroshio  pathway 
was  quite  different.  Compared  to  no  assimilation  or  a  model 
mean,  large  error  reduction  was  obtained  in  all  six  layers  of 
the  model,  including  the  abyssal  layer,  after  about  a  month 
of  simulated  altimeter  data  assimilation.  This  was  tme  even 
when  simulated  data  from  only  one  of  the  altimeters  (any 
one  of  them)  were  assimilated.  The  error  was  significantly 
reduced  by  adding  a  second  altimeter,  with  the  tandem  tracks 
of  Jason- 1  and  TOPEX/Poseidon  clearly  performing  the  best. 
Only  modest  error  reduction  was  obtained  by  adding  a  third 
altimeter  [Hurlburt  et  al. ,  2000;  Smedstad  et  al. ,  2003]. 

2.3 .  Impact  of  Resolution  that  is  Eddy-Resolving  Versus 
Eddy- Permitting  on  Model  Dynamics  and  Simulation  Skill 

Because  of  its  relative  computational  efficiency  and  sim¬ 
plicity,  NTOM  has  been  used  to  perform  a  number  of  studies 
investigating  the  impact  of  increasing  the  horizontal  resolu¬ 
tion  as  well  as  the  impact  of  other  aspects  of  model  design 
on  (a)  model  dynamics,  (b)  simulation  skill,  (c)  convergence 
of  the  mean  circulation  and  variability  in  relation  to  meso- 
scale  features,  and  (d)  agreement  with  observations  in  dif¬ 
ferent  regions  of  the  world  ocean.  Most  of  these  studies  were 
performed  without  ocean  data  assimilation. 

For  the  Kuroshio  near  Japan  (west  of-l  55°E),  Hurlburt  et 
ai  [1996,  1997]  found  that  1/8°  (— 1 5  km)  mid-latitude  reso¬ 
lution  was  sufficient  to  simulate  a  realistic  Kuroshio  pathway, 
including  the  mean  meanders  observed  east  of  Japan.  This 
simulation  included  vertically  compressed  but  otherwise  re¬ 
alistic  bottom  topography.  Corresponding  models  with  1/4° 
resolution  or  with  1/8°  resolution  and  a  flat  bottom  gave  un¬ 
realistic  simulations  of  the  Kuroshio  in  this  region.  The  mean 
Kuroshio  pathway  was  realistic  in  the  1/8°  simulation  with 
topography  because  it  simulated  upper  ocean-topographic 
coupling  via  mixed  baroclinic-barotropic  instabilities  [Hurl¬ 
burt  et  al. ,  1996],  This  coupling  occurs  when  baroclinic 
instability  transfers  energy  to  the  abyssal  layer  where  it  is 
constrained  to  follow  the  geostrophic  contours  of  the  bottom 
topography.  The  eddy-driven  mean  abyssal  currents,  in  turn, 
advect  upper  ocean  current  pathways,  especially  where  these 
currents  intersect  at  large  angles.  In  the  advection  process, 
there  is  often  a  natural  tendency  toward  barotropy  because 
the  current  pathway  advection  is  reduced  when  this  occurs. 
In  the  example  considered  here,  a  modest  seamount  complex 
just  east  of  Japan  played  a  critical  role  in  steering  a  mean 
meander  in  the  Kuroshio  pathway  immediately  east  of  Japan, 
thus  affecting  the  separation  of  the  Kuroshio  from  the  east 
coast  of  Japan. 


Most  of  the  references  in  this  subsection  discuss  upper 
ocean-topographic  coupling  via  flow  instabilities  and  a  the¬ 
ory  for  steering  of  upper  ocean  current  pathways  by  abyssal 
currents,  which  is  discussed  most  comprehensively  by  Hurl¬ 
burt  et  al.  (submitted  manuscript,  2008).  East  of  1 55°E,  Hurl¬ 
burt  et  al  [1996]  and  Hurlburt  and  Metzger  [1998]  found 
that  1/8°  resolution  is  insufficient  and  1/16°  is  sufficient  to 
simulate  the  eastward  penetration  of  the  Kuroshio  as  an  in¬ 
ertial  jet,  the  bifurcation  of  the  Kuroshio  at  the  Shatsky  Rise, 
and  sharp  fronts  that  span  the  northern  Pacific  (the  subarctic 
front  and  a  front  associated  with  the  Kuroshio  Extension). 
Tilburg  et  al.  [2001  ]  also  found  that  1/8°  resolution  is  insuf¬ 
ficient  and  1/16°  is  sufficient  to  simulate  the  observed  mean 
meanders  of  the  Tasman  Front  caused  by  eddy-driven  mean 
abyssal  currents  along  the  slopes  of  ridges  and  trenches  be¬ 
neath  this  front.  Increasing  the  resolution  to  1/32°  gave  only 
modest  improvement. 

Hurlburt  and  Hogan  [2000]  investigated  the  impact  of 
1/8°  (-14  km),  1/16°  (-7  km),  1/32°  (-3.5  km),  and  1/64° 
(~1.8  km)  resolution  on  a  wide  range  of  Gulf  Stream  model 
data  comparisons.  They  found  that  1/16°  resolution  was  the 
minimum  for  realistic  results,  and  substantial  improvement 
was  obtained  by  an  increase  to  1/32°  resolution,  but  only 
modest  improvement  occurred  with  a  further  increase  to 
1/64°  resolution.  These  improvements  include  realistic  Gulf 
Stream  separation  from  the  coast  at  Cape  Hatteras  and  a  real¬ 
istic  pathway  to  the  Grand  Banks  (but  a  more  inertial  current 
at  1/32°  and  1/64°  resolution),  realistic  nonlinear  recirculation 
gyres  that  contribute  to  the  C-shape  of  the  subtropical  gyre, 
realistic  upper  ocean  and  Deep  Western  Boundary  Current 
(DWBC)  transports,  realistic  patterns  and  amplitude  of  SSI  1 
variability  surrounding  the  Gulf  Stream,  realistic  warm  core 
ring  diameters,  population  and  rings  generated  per  year  on  the 
north  side  of  the  Gulf  Stream,  and  realistic  patterns  and  am¬ 
plitude  of  abyssal  EKE.  Hurlburt  and  Hogan  (H.  F.  Hurlburt 
and  P.  J.  Hogan,  The  Gulf  Stream  pathway  and  the  impacts 
of  the  eddy-driven  abyssal  circulation  and  the  Deep  Western 
Boundary  Current,  submitted  to  Dynamics  of  the  Atmosphere 
and  Oceans,  2008]  found  that  steering  by  the  topographically 
constrained,  eddy-driven  mean  abyssal  circulation  alone  is 
sufficient  to  obtain  realistic  Gulf  Stream  separation  from  the 
coast  at  Cape  Hatteras  when  using  1/32°  and  higher  resolu¬ 
tion.  but  assistance  by  the  DWBC  is  required  to  obtain  Gulf 
Stream  separation  in  the  model  with  1/16°  resolution. 

The  Sea  of  Japan  (East  Sea;  JES)  is  characterized  by  weaker 
currents  and  a  smaller  internal  Rossby  radius  of  deformation 
than  the  preceding  regions  [Emery  et  al. ,  1984;  Chelton  et  al. , 
1998;  Oh  et  al. ,  2000).  In  the  JES,  upper  ocean-topographie 
coupling  via  flow  instabilities  has  a  large  impact  on  the 
mean  upper  ocean  circulation,  including  the  pathway  of  the 
subpolar  front  and.  depending  on  the  choice  of  atmospheric 
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forcing  product,  on  the  separation  of  the  East  Korea  Warm 
Current  from  the  Korean  coast.  In  the  resolution  sequence 
1/8°,  1/16°,  1/32°,  1/64°,  this  coupling  does  not  occur  at 
resolutions  coarser  than  1/32°,  but  the  mean  circulation  and 
variability  statistics  are  quite  similar  when  the  resolution  is 
increased  from  1/32°  to  1/64°  [Hogan  and  Hurlburt,  2000, 
2005].  Recently,  anticyclonic  intrathermocline  eddies  (ITEs; 
characterized  by  a  submerged  lens  of  nearly  homogeneous 
water  with  a  bowl-shaped  bottom  and  a  domed  top)  were 
discovered  in  the  JES  [Gordon  et  al .,  2002;  Talley  et  al., 
2004].  Such  eddies  are  far  from  ubiquitous  in  the  world 
ocean  and  were  not  found  in  the  four-layer  NLOM  simu¬ 
lations.  However,  they  were  found  in  15-layer,  1/25°  JES 
HYCOM  simulations  (~3.5-km  resolution,  nearly  the  same 
as  the  1/32°  NLOM;  see  section  2.1),  and  these  simula¬ 
tions  were  used  to  help  explain  why  ITEs  form  in  the  JES 
[Hogan  and  Hurlburt ,  2006].  Like  the  eddies  in  the  1/32° 
NLOM  simulations,  the  eddies  and  ITEs  in  the  JES  tend  to 
be  persistent  and  topographically  constrained  to  a  limited 
range  of  movement.  Further,  mean  SSH,  surface  currents, 
and  abyssal  currents  from  1/25°  JES  HYCOM  and  1/32° 
JES  NLOM  showed  close  quantitative  agreement  when  they 
used  the  same  wind  stress  forcing  [Hurlburt  et  al.,  submitted 
manuscript,  2008]. 

In  a  data-assimilative  example,  Chassignet  et  al.  [2005] 
found  that  the  dynamical  interpolation  skill  of  the  opera¬ 
tional  1/16°  NLOM  prediction  system  was  very  poor  in  the 
northern  Gulf  of  Mexico.  Because  of  insufficient  constraint 
of  the  flow  below  the  sill  depth,  that  system  had  spurious 
southward  leakage  of  a  small  fraction  of  the  DWBC  though 
the  Florida  Strait,  which  then  flowed  cyclonically  around  the 
continental  slope  in  the  Gulf  of  Mexico.  Advection  by  this 
spurious  abyssal  current  is  consistent  with  the  excessively 
rapid  westward  propagation  of  the  upper  ocean  features  in 
the  northern  Gulf  of  Mexico  and  the  poor  dynamical  interpo¬ 
lation  skill  in  the  assimilation  of  altimeter  data  in  that  area. 
The  spurious  DWBC  leakage  was  eliminated  in  the  1/32° 
NLOM  prediction  system,  and  the  dynamical  interpolation 
skill  in  the  northern  Gulf  of  Mexico  was  greatly  improved. 

The  preceding  results  in  this  subsection  and  section  2.2 
point  out  the  striking  differences  between  the  performance 
of  eddy-permitting  and  eddy-resolving  ocean  models  with 
or  without  ocean  data  assimilation.  A  common  theme  is  the 
importance  of  mesoscale  flow  instabilities  in  allowing  eddy- 
driven  abyssal  currents,  constrained  by  bottom  topography, 
to  advect  the  pathways  of  mid-latitude  upper  ocean  currents, 
and  thus  influence  their  mean  pathway.  This  upper  ocean- 
topographic  coupling  requires  that  mesoscale  variability  and 
related  flow  instabilities  be  very  well  resolved  to  obtain  suf¬ 
ficient  coupling,  for  example,  to  achieve  accurate  represen¬ 
tation  of  the  vortex  stretching  and  compression  associated 


with  barochnic  instability  in  the  presence  of  realistic  topog¬ 
raphy.  Thus,  this  eddy-driven  topographic  effect  is  missed  at 
coarser  resolution,  for  example,  in  eddy-permitting  models, 
and  can  lead  to  unexplained  errors  in  simulations  of  mid¬ 
latitude  mean  upper  ocean  current  pathways  and  false  con¬ 
clusions  about  the  influence  of  topography.  Although  based 
on  a  two- layer  theory',  upper  ocean-topographic  coupling  has 
been  demonstrated  in  multilayer  models  using  NL  OM  and 
in  simulations  with  high  vertical  resolution  using  HYCOM 
[Hurlburt  et  al.,  submitted  manuscript,  2008].  This  theory 
applies  in  mid-latitude  regions  with  low  vertical  mode  struc¬ 
ture  in  areas  where  the  topography  does  not  intrude  signifi¬ 
cantly  into  the  stratified  ocean. 

Resolving  the  internal  Rossby  radius  of  deformation  has 
been  used  as  a  criterion  for  distinguishing  between  eddy- 
permitting  and  eddy-resolving  models  because  of  its  relation 
to  the  predominant  space  scale  for  baroclinic  instability.  But 
how  well  does  it  need  to  be  resolved,  and  what  other  crite¬ 
ria  need  to  be  satisfied?  The  preceding  results  provide  some 
very  useful  input.  To  be  eddy-resolving,  a  model  needs  to 
(1)  realistically  simulate  the  eastward  penetration  of  iner¬ 
tial  jets  and  associated  recirculation  gyres  where  they  ex¬ 
ist,  (2)  demonstrate  near  convergence  for  the  mean  strength 
and  pathways  of  both  upper  ocean  and  abyssal  currents  plus 
their  mesoscale  variability  statistics  and,  when  attempting 
to  simulate  the  real  ocean,  it  must  (3)  provide  realistic  com¬ 
parisons  to  observations  of  mean  upper  ocean  and  abyssal 
currents  and  their  variability  statistics,  including  the  basic 
characteristics  of  eddies  and  their  behavior  (adapted  from 
Hurlburt  et  al.,  submitted  manuscript,  2008). 

2.4.  Can  Satellite  Altimetry >  Be  Used  to  Map  Small  Eddies ? 
What  Is  the  Lower  Bound  on  Edd\}  Size  and  its  Dependence 
on  Model  Resolution? 

Plate  2a  shows  a  SeaWiFS  ocean  color  image  in  the  north¬ 
western  Arabian  Sea  and  Gulf  of  Oman.  Such  a  cloud-free 
image  during  a  very  large  outbreak  of  high  chlorophyll  is  an 
unusual  event  in  this  region.  The  image  shows  numerous  ed¬ 
dies  ranging  in  diameter  from  1/4°  to  >  2°,  which  makes  it 
suitable  for  use  in  a  quantitative  evaluation  show  ing  the  ability 
of  satellite  altimeter  data  to  map  mesoscale  variability  when 
assimilated  by  four  global  ocean  models,  models  with  differ¬ 
ing  horizontal  and  vertical  resolution,  model  design,  and  data 
assimilation  techniques,  that  is,  1/16°  NLOM  (8-km  resolu¬ 
tion  at  20°N;  Plate  2b),  1/32°  NLOM  (4  km;  Plate  2c),  1/8° 
NCOM  (18  km;  Plate  2d),  and  1/12°  HYCOM  (8  km;  Plate 
2e).  1/32°  NLOMn  (Plate  2f)  shows  results  from  a  control  run 
designed  to  demonstrate  the  impact  of  the  data  assimilation 
This  run  is  a  repeat  of  the  run  shown  in  Plate  2c,  but  with  no 
assimilation  of  ocean  data.  Based  on  the  results  presented  in 
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Plate  1.  Kuroshio  region  (a,  b)  NAVOCEANO  operational  1/8°  MODAS  SST  analyses  for  (a)  January  L  1999  and  (h)  January 
15,  1999  with  the  color  har  (0.5°C  interval)  designed  to  highlight  the  Kuroshio  pathway,  (c-f )  SSH  analyses  (5-cm  interval)  with 
assimilation  of  altimeter  track  data  from  TOPEX/Poseidon  and  ERS-2  using  the  model  forecast  as  the  lirst  guess  from  (c,  e)  1/16° 
Pacific  NLOM  and  (d,  f)  1/4°  glohal  NLOM  for  (c,  d)  January  1 ,  1999  and  (e,  f)  January  1 5,  1999.  The  altimeter  tracks  assimilated 
on  the  date  given  (using  a  3-day  data  window)  are  overlaid  in  red  on  (c,  e).  (g,  h)  The  14-day  forecast  SSH  from  (g)  1/16°  Pacific 
NLOM  and  (h)  1/4  glohal  NLOM.  Adapted  from  Hurlburt  el  ai  [2000]. 
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Plate  2.  A  comparison  of  eddies  and  currents  seen  in  (a)  ehlorophyll  concentration  from  SeaWiFS  (Oetober  2  6,  2002,  latest  eloud-free 
pixel  composite  with  most  data  from  October  6)  with  (b-e)  data-assimilative  model  nowcast  SSH  and  currents  on  Oetober  6,  2002  from 
(b)  1/16°  global  NLOM  (e)  1/32°  global  NLOM,  (d)  1/8°  global  NCOM,  and  (e)  1/12°  global  HYCOM.  (f)  same  as  (e)  but  without 
ocean  data  assimilation.  In  (a),  eddies  with  clearly-defined  eddy  centers  are  numbered  in  order  of  deereasing  size  for  use  in  Table  2. 
The  color  of  the  number  varies  only  for  visual  elarity  1/32°  NLOM  SSH  and  currents  for  (g)  September  22,  2002;  (h)  September  29, 
2002;  (i)  October  6, 2002  [same  as  (e)]  with  the  most  reeent  week  of  observed  altimeter  tracks  overlaid  as  solid  lines  and  the  remaining 
tracks  as  dashed  lines,  red  for  ERS-2.  black  for  GFO.  and  white  for  Jason- 1.  At  the  time.  Jason- 1  data  were  not  available  in  real  time  tit 
NAVOCEANO  and  thus  not  assimilated  by  the  operational  1/16°  NLOM  system  which  provided  the  results  used  in  (b).  Note  that  a  few 
of  the  smallest  eddies  are  not  visible  due  to  limitations  of  the  small  plots  or  beeause  they  were  slightly  east  of  63°E  in  a  model  result. 
NLOM  results  are  adapted  from  Shriver  et  al.  [2007]. 
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section  2.3,  with  18-km  resolution  at  20°N  and  ~15  km  at 
mid-latitudes,  the  1/8°  NCOM  system  is  eddy-permitting,  and 
the  rest  of  the  systems  are  eddy-resolving.  All  of  the  model 
runs  include  atmospheric  forcing;  all  of  the  data-assimilative 
runs  include  assimilation  of  SST,  and  1/12°  HYCOM  also 
assimilates  T  &  S  profiles,  but  satellite  altimetry  is  the  key 
observ  ing  system  for  mapping  the  eddies.  Results  of  this  com¬ 
parison  for  1/16°  and  1/32°  NLOM  can  be  found  in  Shriver 
et  al.  [2007].  All  of  the  systems  assimilated  altimeter  data 
from  the  ERS-2,  GFO,  and  Jason- 1  missions,  except  that  the 
operational  1/16°  NLOM  system  omitted  Jason- 1  because  the 
data  were  not  yet  available  in  real  time  at  NAVOCEANO. 
Shriver  et  al.  [2007]  found  that  the  Jason- 1  data  added  only 
minimal  value  in  this  particular  case,  which  is  consistent  with 
the  results  of  Smedstad et  al.  [2003],  who  used  simulated  data 
and  found  that  adding  a  third  altimeter  provided  only  modest 
error  reduction. 

A  quantitative  analysis  of  model  eddy  center  location 
error  in  comparison  to  the  ocean  color  image  (Plate  2a)  is 
shown  in  Table  2.  All  eddies  with  a  clearly  defined  eenter 
were  used  except  for  a  cyclonic  eddy  located  inside  anticy- 
elonie  eddy  3.  The  eddy  eenter  locations  in  the  models  and 
the  oeean  color  were  determined  as  independently  as  pos¬ 
sible  and  compared  only  to  the  extent  needed  to  match  them 
up  (location,  rotation  direction,  and  size).  The  accuracy  in 
determining  the  locations  of  the  eddy  centers  is  typically 
IO  I5  km  for  both  the  ocean  color  and  the  models  and  thus 
is  15  20  km  in  determining  the  eddy  center  location  error. 
The  eddies  in  Table  2  are  listed  in  the  order  of  decreasing 
diameter,  as  seen  in  the  ocean  color,  to  assist  in  determining 
the  minimum  eddy  size  that  eould  be  mapped  by  assimila¬ 
tion  of  the  altimeter  data  in  each  model.  The  first  four  eddies 
are  ~2°  or  more  in  diameter,  eddies  5  8  are  ~1°,  12  16  are 
~l/2°,  and  17  20  are  -1/4°.  With  a  few  exceptions,  model 
eddies  1 1  through  20  are  typically  ~50%  or  more  larger  than 
seen  in  the  ocean  color,  but  the  rest  are  generally  consist¬ 
ent  with  the  oeean  eolor  imagery  Using  model-independent 
analyses  of  TOPEX/Poseidon  and  ERS-I  and  2  altimeter 
data,  Ducet  et  al.  [2000]  estimated  - 75  km  as  the  minimum 
eddy  diameter  that  could  be  mapped  using  satellite  altimeter 
data.  Thus,  75  km  was  used  to  divide  the  observed  eddies 
into  large  and  small  eddies,  making  eddies  1 1  through  20  the 
small  eddies. 

Table  2  summarizes  the  results  of  comparing  eddy  center 
locations  in  the  models  with  those  in  the  ocean  color.  In  this 
comparison,  the  ability  to  map  eddies  as  small  as  1/4°  in 
diameter  is  not  limited  by  the  satellite  altimetry,  but  by  the 
model  resolution.  The  1/32°  NLOM  system  depicted  90% 
of  the  eddies,  80%  of  the  small  eddies,  and  all  four  of  the 
1/4°  eddies  (three  of  the  1/4°  eddies  with  an  accurate  loca¬ 
tion)  versus  the  1/8°  NCOM  system,  which  depieted  55%  of 


the  eddies,  80%  of  the  larger  eddies  (I  10)  and  30%  of  the 
small  eddies.  Note  that  the  1/8°  NCOM  system  assimilated 
the  sterie  SSI  I  anomalies  from  the  1/32°  NLOM  system  in 
the  form  of  synthetic  T  &  S  profiles  [Rhodes  et  al .,  2002 
and  section  1.2].  Generally,  about  seven  grid  intervals  was 
the  minimum  eddy  diameter  the  models  could  realistically 
depiet,  that  is,  128  km  for  1/8°  NCOM,  56  km  for  1/12°  HY- 
COM  and  1/16°  NLOM,  and  28  km  for  1/32°  NLOM.  In 
some  cases,  eddies  were  depicted  by  increasing  the  diameter 
of  small  eddies  to  at  least  this  minimum  size,  for  example, 
eddy  11  (but  not  eddy  18)  in  1/8°  NCOM  and  eddy  18  in 
1/16°  NLOM.  Consistent  with  the  seven  grid  interval  mini¬ 
mum,  1/16°  NLOM  was  able  to  map  eddies  1/2°  and  larger, 
but  not  the  1/4°  eddies.  All  of  the  eddies  seen  in  the  ocean  eolor 
were  mapped  by  at  least  one  of  the  data-assim dative  systems. 

In  comparison  to  I  / 32°  NLOMn  (with  no  ocean  data  assim¬ 
ilation),  the  impact  of  the  altimeter  data  is  strikingly  evident 
(Plate  2  and  Table  2).  As  1/32°  NLOMn  depicts  numerous 
small  eddies  (Plate  20  and  five  of  these  match  up  with  50% 
of  the  small  eddies  seen  in  the  ocean  eolor,  one  might  ques¬ 
tion  the  impact  of  the  altimeter  data  assimilation  in  mapping 
these  eddies.  However,  Table  2  clearly  shows  the  impact  of 
the  assimilation  on  the  depiction  of  small  eddies,  for  exam¬ 
ple,  by  comparing  1/32°  NLOM  and  1/32°  NLOMn,  we  find 
that  the  depiction  of  small  eddies  is  80%  versus  50%,  the 
more  accurate  eddy  center  location  is  67%  versus  33%,  the 
median  eddy  center  location  error  is  22.5  versus  42  km,  and 
the  minimum  eddy  center  location  error  is  I  1  versus  26  km. 
Furthermore,  62.5%  (87.5%)  of  the  small  eddies  in  1/32° 
NLOM  have  eddy  center  location  errors  less  than  the  mini¬ 
mum  of  26  km  (median  of  42  km)  in  1/32°  NLOMn 

Clearly,  the  minimum  eddy  diameter  mapped  using  satel¬ 
lite  altimetry  is  much  less  than  75  km  when  the  altimeter 
data  are  assimilated  by  an  cddy-rcsolving  ocean  model  (all 
except  the  eddy-pernutting  1/8°  NCOM  in  Plate  2  and  Ta¬ 
ble  2).  The  1/32°  NLOM  system  (Plates  2g  2i)  is  used  to 
illustrate  how  a  model  can  assist  in  the  depiction  of  eddies 
and  currents  by  using  its  dynamical  interpolation  skill.  The 
altimeter  tracks  crossing  the  region  are  overlaid  on  each 
panel  (the  most  recent  week  as  solid  lines  and  the  rest  as 
dashed  lines),  although  only  a  3-day  data  window  is  used 
in  the  daily  assimilation  cycle.  In  total,  these  panels  show 
the  most  recent  three  weeks  of  altimeter  data,  September  16 
through  October  6,  2002. 

A  striking  example  of  model  dynamical  interpolation  skill 
is  provided  by  a  current  that  enters  the  southern  boundary  of 
the  region  near  60°E  and  ultimately  wraps  around  eddy  12. 
It  is  marked  by  a  ribbon  of  blue  (lower  chlorophyll)  in  the 
oeean  color  and  flows  northward  along  -60°  E,  then  west¬ 
ward  north  of  1 8°N  and  then  northward  along  ~58.5°E  until  it 
wraps  around  eddy  12,  which  is  1/2°  in  diameter.  This  current 
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is  seen  in  all  of  the  data-assimilative  models.  Eddy  12  and 
adjacent  eddy  1 3  are  depicted  in  all  the  data-assimilative  mod¬ 
els  except  for  eddy  13  in  1/8°  NCOM.  The  non-assimilative 
1/32°  NLOMn  depicts  neither  the  current  nor  the  eddies.  Dur¬ 
ing  the  week  leading  up  to  October  6,  2002,  the  entire  system 


is  unobserved  by  the  altimetry  except  for  the  southernmost 
portion  of  the  current.  The  system  is  quite  well  observed  dur¬ 
ing  September  \(y  22,  when  it  is  significantly  different.  From 
September  16  to  October  6,  eddies  12  and  13  are  never  more 
than  peripherally  observed,  and  eddy  13  does  not  even  exist 


Table  2.  Eddy  Center  Location  Errors  in  Ocean  Prediction  Models  Compared  to  Ocean  Color  From  SeaWiFS  in  the  Northwestern  Arabian 
Sea  and  Gulf  of  Oman. 


Ocean  Color 
Eddy  ID# 

A  or  C 

Ocean  Color  Eddy 
Center  Location 

1/16°  NLOM 

1/32°  NLOM 

1/8°  NCOM 

1/12°  HYCOM 

1/32°  NLOMn 
No  Assim 

°N 

°E 

Eddy  Center  Position  Error,  km 

1 

C 

23.55 

60.2 

35 

18 

80 

48 

NP 

2 

A 

18.3 

62.0 

103 

28 

66 

48 

NP 

3 

A 

21.15 

60.25 

44 

58 

48 

37 

NP 

4 

C 

19.65 

60.4 

43 

12 

15 

59 

45 

5 

C 

16.8 

55.65 

35 

42 

75 

72 

31 

6 

c 

17.0 

58.9 

42 

17 

38 

68 

NP 

7 

c 

16.7 

57.9 

53 

79 

NP 

97 

NP 

8 

A 

18.0 

57.6 

NP 

40 

36 

68 

NP 

9 

c 

25.1 

57.6 

NP 

39 

91 

76 

NP 

10 

A 

24.7 

62.5 

30 

35 

NP 

NP 

NP 

1  1 

C 

23.7 

62.3 

a 

22 

37 

52 

42 

12 

c 

19.3 

58.8 

30 

30 

50 

65 

NP 

13 

A 

19.225 

59.35 

35 

11 

NP 

26 

NP 

14 

A 

25.3 

58.55 

36 

33 

NP 

28 

30 

15 

C 

24.1 

61.75 

55 

NP 

NP 

NP 

47 

16 

C 

22.25 

62.7 

14 

* 

NP 

48 

- 

17 

c 

23.1 

62.55 

NP 

13 

NP 

44 

NP 

18 

A 

22.5 

62.85 

18 

51 

47 

NP 

44 

19 

c 

22.05 

62.05 

NP 

23 

NP 

NP 

26 

20 

c 

22.2 

59.95 

NP 

12 

NP 

22 

NP 

%  of  Eddies  Present  in  the  Model 

All  eddies 

70 

90 

55 

80 

35 

Large  eddies,  1 

1-10 

80 

100 

80 

90 

20 

Small  eddies,  1 

11-20 

60 

80 

30 

70 

50 

Median  Eddy  Center  Position  Error,  km 

All  eddies 

35.5 

29 

48 

50 

42 

Large  eddies 

42.5 

37 

57 

68 

38 

Small  eddies 

32.5 

22.5 

47 

44 

42 

%  of  eddies  With  Most  Accurate  Position 

All  eddies 

22.5 

52.5 

5 

10 

10 

Large  eddies 

20 

50 

10 

10 

10 

Small  eddies 

25 

55 

0 

10 

10 
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on  September  22  and  29.  A  GFO  track  does  appear  to  play 
an  important  role  in  narrowing  the  westward  bulge  of  eddy 
2  from  September  22  to  29,  which  assists  in  the  detachment 
of  eddy  13.  However,  model  dynamics  are  required  to  move 
eddy  13  northward  and  partially  merge  it  with  an  existing 
anticyclonic  eddy  north  of  eddy  12.  This  process  gave  an  I  1- 
km  eddy  center  location  error  for  eddy  13  in  the  1/32°  NLOM 
system. 

In  a  second  example,  eddy  19  was  well  observed  by  an 
ERS-2  track  during  the  week  leading  up  to  September  22. 
The  result  was  a  strong  small  eddy  in  1/32°  NLOM  and  a 
weak  one  in  1/16°  NLOM.  The  1/32°  NLOM  system  main¬ 
tained  this  nearly  stationary  eddy  through  October  6,  but  it 
quickly  dissipated  in  the  1/16°  NLOM  system.  During  the 
week  leading  to  October  6,  eddy  16  was  well  observed  by 
GFO  and  ERS-2  altimeter  tracks,  and  it  was  represented  in 
1/16°  NLOM  with  a  14-km  eddy  center  location  error.  Eddy 
16  was  not  represented  in  1/32°  NLOM  because  eddy  19 
was  too  large  and  the  observations  of  eddy  16  were  treated 
as  observations  of  the  eastern  edge  of  eddy  19,  thereby  mak¬ 
ing  it  a  fusion  of  eddies  16  and  19. 

Cyclonic  eddy  20  is  the  smallest  of  the  1/4°  eddies  (24  km 
in  diameter  in  the  1/32°  NLOM  system).  It  lies  adjacent  to 
the  coast  and  a  strong  current  that  separates  from  the  coast. 
In  the  ocean  color  image,  this  current  wraps  around  large 
anticyclonic  eddy  3  Eddy  20  was  not  observed  by  altimetry, 
but  forms  in  1/32°  NLOM  and  1/12°  HYCOM  due  to  the 
current  separation  from  the  boundary  and  proximity  to  the 
boundary.  This  small  eddy  does  not  form  in  the  coarser  reso¬ 
lution  1/16°  NLOM  and  1/8°  NCOM,  although  they  depict 
eddy  3  and  the  separating  boundary  current  with  comparable 
realism.  The  separating  boundary  current  is  shown  with  less 


accuracy  in  1/12°  HYCOM,  but  HYCOM  does  simulate  a 
cyclonic  eddy  similar  to  eddy  20  adjacent  to  the  coast. 

The  preceding  results  demonstrate  that  an  eddy-resolving 
ocean  model  can  be  skillful  in  assimilating  data  along  multiple 
altimeter  tracks  from  multiple  satellites  and  in  dynamically  in¬ 
terpolating  it  through  time  to  build  coherent  mesoscale  features 
that  are  consistent  with  observations,  including  small  eddies. 
They  also  demonstrate  that  increasing  the  horizontal  resolution 
of  the  model  increases  that  skill,  especially  in  mapping  small 
eddies  25  75  km  in  diameter  In  addition,  the  results  show  that 
a  model  with  dynamical  interpolation  skill  can  use  the  data 
assimilation  to  produce  an  environment  more  conducive  to 
the  formation  of  poorly  observed  eddies  than  the  same  model 
without  the  assimilation,  a  capability  that  is  especially  useful  in 
representing  small  eddies,  as  demonstrated  here. 

2.5.  Eddy- Resolving  Ocean  Model  Forecast  Skill  and  the 
Impacts  of  Atmospheric  Forcing ,  Model  Choice ,  and  Model 
Resolution 

Several  metrics  are  used  routinely  to  assess  model  fore¬ 
cast  skill  based  on  comparison  to  a  verifying  analysis.  To 
demonstrate  useful  skill,  a  forecast  must  be  superior  to  per¬ 
sistence  (a  forecast  of  no  change)  and  to  either  climatology 
or  some  minimum  value  of  the  metric.  The  metrics  are  root 
mean  square  (rms)  error,  anomaly  correlation  (AC),  skill 
score,  and  for  the  Kuroshio  and  the  Gulf  Stream,  axis  error. 
See  Smedstad  et  ai  [2003]  for  discussion  and  examples  of 
results  from  each  of  these.  Normalized  mis  error,  normalized 
by  the  standard  deviation  of  the  anomalies,  is  another  useful 
forecast  metric.  Here,  only  the  SSH  is  evaluated  and  only  the 
AC  is  used  a  metric  commonly  used  in  verify  ing  weather 


Notes  to  Table  2: 

aA  single  “fused”  model  eddy  represents  two  eddies  in  the  ocean  color  and  the  model  eddy  center  lies  between  the  two  observed  eddies. 
The  model  eddy  center  position  error  is  listed  only  under  the  closer  ocean  color  eddy. 

A,  anticyclonic  eddies;  C,  cyclonic  eddies;  NP,  eddy  not  present 

The  ocean  color  eddy  11)  numbers  are  plotted  in  Plate  2a.  Eddies  are  listed  in  order  of  decreasing  size  as  depicted  by  the  ocean  color  Eddy 
position  measurement  error  is  10-15  km  in  both  the  ocean  color  and  the  models. 

1/32°  NLOM;  1/32°  global  NLOM  7-layer  prediction  system  with  assimilation  of  altimeter  track  data  from  ERS-2,  GEO  and  Jason- 1 
altimeters  (presently  operational  at  NAVOCEANO). 

1/16°  NLOM;  Then  operational  1/16°  global  NLOM  seven  layer  prediction  system  with  assimilation  of  real-time  altimeter  track  data  from 
the  ERS-2  and  GFO  altimeters  (Jason- 1  was  not  in  the  operational  data  stream  at  that  time) 

1/8°  NCOM;  1/8°  global  NCOM  40-level  prediction  system  with  assimilation  of  l/32‘  NLOM  SSH  via  synthetic  T  &  S  profiles  (presently 
operational  at  NAVOCEANO). 

1/12°  HYCOM.  .08°  global  HYCOM  prediction  system  with  32  hybrid  layers  and  assimilation  of  altimeter  track  data  from  ERS-2,  GFO 
and  Jason- 1  altimeters  plus  T  &  S  profiles  (presently  pre-opcrational  and  running  in  real  lime  at  NAVOCEANO). 

1/32°  NLOMn:  1/32°  global  NLOM  with  no  assimilation  of  ocean  data,  only  atmospheric  forcing. 

Some  model  versions  of  eddy  1 8  lie  east  of  63°E  and  thus  are  not  depicted  in  Plate  2. 
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forecasts.  It  is  calculated  with  the  same  methodology  used  by 
Smedstad  et  ai  [2003]  and  Shriver  et  ai  [2007]: 

YV-fXa-a) 

AC(fa)=  - - 

where  /is  the  forecast  value  and  a  is  the  analysis  value  at 
a  given  time  and  model  grid  point.  Here,  the  mean  of  the 
forecast /  and  the  analysis  a  are  identical,  that  is,  the  modi¬ 
fied  model  mean  SSH  used  in  assimilating  SSH  anomalies 
observed  by  satellite  altimetry. 

Some  of  the  real-time  forecast  models  perform  30-day 
forecasts  once  a  week.  For  these  forecasts,  the  real-time  at¬ 
mospheric  forcing  reverts  toward  climatology  after  the  end 
of  the  forecast  atmospheric  forcing  (after  5  days  in  these  ex¬ 
amples).  To  evaluate  the  impact  of  real-time  versus  analysis- 
quality  atmospheric  forcing  on  ocean  model  forecast  skill, 
two  sets  of  “forecasts”  were  performed  and  evaluated  using 
1/12°  global  HYCOM  and  1/32°  global  NLOM,  but  only 
real-time  forcing  in  the  1/16°  global  NLOM  forecasts.  In  ad¬ 
dition,  the  SSH  forecasts  were  verified  against  tide  gauge 
data,  an  independent,  unassimilated  data  set. 

The  results  of  forecast  verification  versus  analyses  and 
tide  gauge  data  are  presented  in  Plate  3.  The  global  forecast 
verification  covers  only  regions  with  assimilation  of  altim¬ 
eter  data,  from  50°S  to  65°N  for  NLOM  and  603S  to  47°N  in 
these  early  results  for  HYCOM.  The  Gulf  Stream  subregion 
was  chosen  to  illustrate  forecast  skill  in  a  highly  nonlinear 
region  with  strong  flow  instabilities,  current  meanders,  and 
generation  of  mesoscale  eddies.  These  are  class  2  (Table 
1)  phenomena  which  are  insensitive  to  errors  in  the  atmos¬ 
pheric  forcing  on  a  30-day  timescale.  T  he  equatorial  Pacific 
has  significant  elements  from  all  four  classes,  including  class 
1  phenomena  such  as  the  onset  of  wind-driven  equatorial 
waves  and  variations  in  equatorial  upwelling.  Thus,  30-day 
forecasts  in  this  region  are  more  sensitive  to  errors  in  the 


atmospheric  forcing.  The  region  in  the  NW  Arabian  Sea  and 
Gulf  of  Oman  was  chosen  to  show  forecast  skill,  a  prerequi¬ 
site  for  dynamical  interpolation  skill,  over  nearly  the  same 
region  shown  in  Plate  2.  The  forecast  results  are  quite  typical 
of  the  world  ocean  with  modest  sensitivity  to  the  errors  in 
atmospheric  forcing  beyond  the  first  2  weeks. 

A  dramatically  different  result  is  seen  in  the  forecast  re¬ 
sults  for  the  Persian  Gulf  from  the  1/12°  global  HYCOM 
system.  Here,  persistence  demonstrates  forecast  skill  (AC 
>0.6)  for  only  2  days  and,  similarly,  the  forecast  with  real¬ 
time  forcing  lasts  only  2  days  beyond  the  end  of  the  2  days 
of  forecast  forcing  (represented  by  5  days  of  analysis-quality 
forcing  in  these  tests).  In  contrast,  forecast  skill  remains  high 
for  the  30-day  duration  of  the  forecasts  when  analysis  qual¬ 
ity  forcing  is  used.  The  Persian  Gulf  is  a  shallow-water  basin 
where  the  SSH  exhibits  a  largely  deterministic  class  1  response 
to  atmospheric  forcing.  As  a  result,  the  forecast  skill  is  char¬ 
acterized  by  rapid  loss  of  initial  state  impact,  and  it  is  strongly 
determined  by  the  timescale  of  accurate  atmospheric  forcing 
(i.e.,  it  is  limited  by  the  timescale  of  atmospheric  predictive 
skill  in  a  real-time  prediction  system).  Similar  results  were  ob¬ 
tained  for  the  shallow  Yellow  and  Bohai  Seas  region  north  of 
30°N.  Shriver  et  al.  [2007]  show  forecast  results  from  more 
subregions,  using  the  1/32°  global  NLOM  system,  including 
results  for  two  deep-water  semi-enclosed  seas.  For  these,  the 
forecast  skill  is  more  characteristic  of  other  deep-water  regions. 

Overall,  the  30-day  forecast  results  show  low  sensitivity 
to  the  difference  between  the  climatological  and  analysis- 
quality  atmospheric  forcing  fields  in  deep-water  regions. 
Outside  the  surface  mixed  layer  and  shallow  water,  the  evo¬ 
lution  of  the  forecasts  is  more  sensitive  to  the  initial  state 
than  to  the  atmospheric  forcing  anomalies  over  most  of 
the  global  ocean.  Much  of  this  is  due  to  the  fact  that  meso¬ 
scale  eddies  are  not  confined  to  highly  energetic  regions  of 
the  world  ocean  like  the  Gulf  Stream  and  Kuroshio,  but  are 
nearly  ubiquitous  in  the  world  ocean,  as  demonstrated  by 
model-independent  analyses  of  satellite  altimeter  data  [Le 


Plate  3.  (Opposite)  Verification  of  30-day  ocean  forecasts,  (a-i)  median  SSH  anomaly  correlation  versus  forecast  length  in  comparison 
with  the  verifying  analysis  for  (a,  e),  the  global  domain  (45°S-45°N)  and  four  subregions  (b-d,  f-i)  defined  in  Tabic  3.  (a-d)  are  for  1/32° 
NLOM  (solid  lines)  and  1/16°  NLOM  (dashed  lines)  and  (e-i)  are  for  1/12°  HYCOM  The  red  curves  verify  forecasts  using  operational 
atmospheric  forcing  which  reverts  toward  climatology  after  five  days.  The  green  curves  verify  “forecasts”  with  analysis-quality  forcing  for 
the  duration,  and  the  blue  curves  verify  forecasts  of  persistence  (i.e.,  no  change  from  the  initial  stale).  The  plots  for  NLOM  give  the  median 
statistics  over  twenty-two  30-day  forecasts  initialized  during  the  period  from  June  1 , 2001  through  May  3 1 , 2002,  and  for  HYCOM  twenty 
30-day  forecasts  initialized  during  the  period  from  January  2004  through  December  2005,  both  periods  when  data  from  three  nadir-beam 
altimeters  in  corresponding  orbits  were  assimilated.  The  same  forecast  results  were  used  to  obtain  (j-1)  median  correlation  between  forecast 
and  observed  SSH  fluctuations  from  1/12°  HYCOM  with  operational  forcing  during  the  forecast  (red  lines),  1/12°  HYCOM  with  analysis 
quality  forcing  for  the  duration  (green  lines),  1/16°  (blue  lines),  and  1/32°  (black  lines)  NLOM  (both  with  operational  atmospheric  forcing) 
at  (j)  23  (49)  open  ocean  island  tide  gauge  stations  for  HYCOM  (NLOM),  (k)  91  (29)  coastal  tide  gauges  for  HYCOM  (NLOM),  and  (I)  all 
1  14  (78)  tide  gauges  for  HYCOM  (NLOM).  A  1 3-day  moving  average  was  applied  to  filter  time  scales  not  resolved  by  the  altimeter  data. 
Tide  gauge  SSH  data  are  not  assimilated  by  the  ocean  prediction  systems.  NLOM  results  are  adapted  from  Shriver  et  al.  (2007). 
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Table  3.  Subregions  Depicted  in  Plate  3 


Subregion 

Latitude  Range 

Longitude  Range 

World  Ocean 

45°S-45°N 

All 

Gulf  Stream  region 

35°N-45°N 

76°W-40°W 

Equatorial  Pacific 

20°S-20°N 

I09°E-77°W 

NW  Arabian  Sea 
and  Gulf  of  Oman 

15°N-26°N 

5 1 °E-65°E 

Traon  et  al. ,  1998]  and  as  illustrated  in  Plate  2a  for  a  range 
of  space  scales  in  a  region  without  high  SSH  variability. 

Comparing  AC  versus  time  from  the  30-day  forecasts  per¬ 
formed  using  1/16°  NLOM  and  1/32°  NLOM,  one  finds  the 
perhaps  surprising  result  that  globally  and  for  many  subregions, 
the  AC  is  higher  in  the  1/16°  model  [Shriver  et  al. ,  2007],  The 
higher  resolution  model  has  more  smaller  scale  features  with 
larger  amplitude  that  can  get  out  of  phase  on  shorter  times¬ 
cales.  Consistent  with  this  difference  between  the  two  models, 
persistence  also  has  less  skill  in  the  1/32°  NLOM  system  than 
in  the  1/16°  system.  Thus,  the  best  way  to  judge  performance 
is  not  by  the  AC  for  the  forecast,  but  by  the  spread  between 
the  AC  of  the  forecast  and  persistence.  For  a  given  day,  that  is 
always  greater  for  the  1/32°  NLOM  system. 

The  preceding  results  also  highlight  the  need  to  assess 
forecast  skill  using  comparisons  to  independent,  unassimi¬ 
lated  data  sets.  Thus,  Plates  3j  31  show  evaluation  of  SSH 
from  1/12°  global  HYCOM  and  the  two  global  NLOM  sys¬ 
tems  by  comparison  to  unassimilated  tide  gauge  data.  For 
HYCOM  that  also  includes  a  comparison  of  “ forecast s”  us¬ 
ing  operational  and  analysis-quality  forcing.  A  13-day  mov¬ 
ing  average  was  applied  to  filter  timescales  not  resolved  by 
the  assimilated  altimeter  data.  Again,  these  comparisons 
indicate  that  the  1/32°  NLOM  system  has  greater  forecast 
skill  than  the  1/16°.  HYCOM  forecasts  give  greater  skill 
at  coastal  stations  (Plate  3k)  than  island  stations  (Plate  3j), 
while  the  reverse  is  true  for  the  NLOM  systems  HYCOM 
forecasts  demonstrate  greater  sensitivity  to  analysis-quality 
versus  operational-quality  atmospheric  forcing  at  coastal 
stations  than  at  island  stations,  as  expected  from  the  differ¬ 


ing  classes  of  response  to  atmospheric  forcing  in  coastal  ver¬ 
sus  deep  ocean  regions  (Table  2).  HYCOM  outperformed 
NLOM  at  coastal  stations  because  it  includes  shallow  water 
and  because  the  HYCOM  numerical  grid  generally  extends 
closer  to  the  tide-gauge  locations.  In  contrast,  NLOM  out¬ 
performed  HYCOM  at  island  stations  on  forecast  timescales 
longer  than  2  weeks  where  model  grid  proximity  to  the  tide 
gauges  is  similar  for  both  models.  In  addition,  HYCOM 
forecast  accuracy  decreased  more  rapidly  than  NLOM’s  at 
island  stations.  As  the  HYCOM  and  NLOM  forecasts  were 
performed  for  different  years  and  there  are  large  differences 
in  the  tide  gauge  stations  used,  only  large  HYCOM-NLOM 
differences  were  considered  in  the  preceding  comparisons 
between  the  systems. 

2.6 .  Mean  SSH  for  Assimilation  of  Altimeter  Data:  Roles  of 
Model  and  Observation  Based  Means 

The  capability  to  simulate  mean  SSH  accurately  is  impor¬ 
tant  for  an  ocean  model  because  (1)  an  accurate  mean  SSH 
field  is  required  for  addition  to  the  deviations  provided  by 
satellite  altimetry,  and  a  sufficiently  accurate  model  can  be 
used  in  providing  it,  and  (2)  if  the  mean  SSH  is  not  well 
simulated  by  the  model,  the  model  forecasts  will  tend  to  drift 
with  a  bias  toward  the  erroneous  model  mean,  reducing  the 
AC  scores  of  the  forecasts  (discussed  in  section  2.5).  In  a 
data-assimilative  ocean  model,  high  accuracy  in  SSH  fields 
is  essential  to  properly  represent  ocean  currents  and  fronts 
and  to  avoid  contributing  to  errors  in  subsurface  T  &  S. 
Thus,  accurate  mean  SSH  is  a  critical  application  of  a  class  3 
ocean  response  to  atmospheric  forcing  in  ocean  data  assimi¬ 
lation  Model  and  observation-based  means  each  have  their 
advantages  and  disadvantages  for  this  purpose.  Fortunately, 
the  two  types  of  mean  are  converging.  The  standard  devia¬ 
tion  of  the  difference  between  the  observation-based  Maxi¬ 
menko  and  Niiler  [2005]  SSH  mean  (Plate  4a)  and  a  1/12° 
global  HYCOM  mean  (Plate  4b)  is  only  9.3  cm,  and  the  bias 
between  their  global  areal-averaged  means  is  7.1  cm.  The 
detailed  agreement  of  the  frontal  segments  between  30°S 
and  60°S  is  particularly  striking. 


Plate  4.  (Opposite)  Mean  SSH  (a)  from  Maximenko  and  Niiler  [2005]  based  on  satellite  altimeter,  drifting  buoy  and  wind  data  spanning 
1992  2002  and  an  improved  geoid,  GRACE  (Gravity  Recovery  and  Climate  Experiment)  Gravity  Model  01  (GGM01)  and  (b)  simulated 
by  1/12°  global  HYCOM  using  ECMWF  climatological  wind  and  thermal  forcing.  The  standard  deviation  of  the  difference  between  them 
is  9.3  cm,  and  the  bias  between  their  global  areal -averaged  mean  is  7.1  cm.  (e)  MOD  AS  mean  sterie  height  anomaly  relative  to  1 ,000  m 
based  solely  on  historical  hydrography.  Zoom  views  of  (d)  the  Maximenko  and  Niiler  mean  SSH  and  (c)  the  I /1 2°  global  HYCOM  mean 
SSH  in  the  western  equatorial  Pacific  and  adjacent  marginal  seas,  (f,  g)  Mean  currents  in  the  northeast  Pacifie  with  mean  current  speed  in 
color  from  1/8°  global  NCOM  (f)  w  ithout  ocean  data  assimilation  and  (g)  with  assimilation  of  SSH  anomalies  analyzed  by  NLOM  from 
satellite  altimetry  using  the  MODAS  mean  SSI  I  in  (e).  Without  the  assimilation,  the  Alaska  Stream,  marked  b\  the  ribbon  of  red  in  (f), 
is  robustly  simulated  by  1/8°  global  NCOM,  but  it  is  greatly  weakened  in  (g)  when  the  data  assimilation  is  added  because  the  mean  SSH 
associated  with  the  Alaska  Stream  is  poorly  represented  in  the  MODAS  mean.  Panels  (f)  and  (g)  are  adapted  from  Barron  et  al  [2007). 
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The  HYCOM  simulation  was  initialized  from  a  hydro- 
graphic  climatology,  and  the  mean  over  years  9-13  is  shown 
in  Plate  4b.  Although  ocean  data  assimilation  could  have 
been  included,  this  simulation  used  only  climatological  at¬ 
mospheric  forcing,  except  for  weak  relaxation  to  sea  surface 
salinity  (in  addition  to  surface  salinity  forcing  from  evapo¬ 
ration,  precipitation,  and  rivers).  The  monthly  atmospheric 
climatology  is  based  on  the  European  Centre  for  Medium- 
Range  Weather  Forecasts  (ECMWF)  reanalysis  (ERA- 1 5) 
[Gibson  et  al .,  1 997]  with  the  1 0-m  winds  converted  to  wind 
stress  using  the  bulk  formulation  of  Kara  et  al.  [2005].  Sub¬ 
monthly  fluctuations  were  added  to  the  wind  stress  to  ener¬ 
gize  the  mixed  layer. 

The  MODAS  mean  steric  height  anomaly  relative  to 
1 ,000  m  (Plate  4c)  is  based  solely  on  historical  hydrography. 
Such  means  have  insufficient  data  to  accurately  map  sharp 
features  in  the  mean,  such  as  boundary  currents  and  ocean 
fronts.  This  deficiency  causes  serious  problems  in  ocean  data 
assimilation  as  illustrated  in  Plates  4f  and  4g  showing  mean 
currents  in  the  NE  Pacific  from  1/8°  global  NCOM.  Plate  4f 
shows  that  without  ocean  data  assimilation,  NCOM  simulates 
a  robust  Alaska  Stream  (red  ribbon  of  high  current  speed), 
a  current  with  realistic  transport  in  comparison  to  9  years 
of  observations  presented  by  Onishi  and  Ohtani  [1999]  and 
Onishi  [2001].  However,  when  the  MODAS  mean  is  used 
as  the  reference  tor  assimilation  of  altimetric  anomalies,  the 
Alaska  Stream  in  NCOM  is  greatly  weakened.  This  problem 
was  rectified  by  using  a  modified  model  mean. 

Maximenko  and  Niiler  [2005]  derived  a  1/2°  absolute  SSH 
mean  (i.e.,  not  relative  to  a  reference  depth)  based  on  data  from 
drifting  buoys,  satellite  altimetry,  winds,  and  the  GRACE 
(Gravity  Recovery  and  Climate  Experiment)  Gravity  Model 
01  (GGM01)  geoid  [see  also  Niiler  et  al .,  2003].  It  tends  to 
have  more  sharply  defined  features  than  the  MODAS  mean 
surface  dynamic  height  relative  to  1000  m  from  hydrography. 
However,  mean  fronts  and  boundary  currents  tend  to  be  the 
most  sharply  defined  in  the  model  mean.  In  addition,  some 
boundary  currents  are  poorly  represented  in  the  Maximenko 
and  Niiler  mean,  for  example,  the  Mindanao  Current  along  the 
east  coast  of  the  southern  Philippines  (Plate  4d  compared  to 
the  model  mean  shown  in  Plate  4e)  and  the  unrealistic  features 
around  the  island  of  Borneo. 

Thus,  it  is  desirable  to  use  a  mean  SSH  field  from  an  ac¬ 
curate  eddy-resolving  ocean  model  that  was  initialized  from 
climatology.  Mean  SSH  features  have  a  wide  range  of  am¬ 
plitudes  and  space  scales,  and  an  eddy-resolving  model 
should  have  an  advantage  in  representing  mean  SSH  associ¬ 
ated  with  smaller  gyres  and  weaker  currents  and  features  as¬ 
sociated  w  ith  complex  geometry,  such  as  the  Indonesian  and 
Philippine  archipelagos  (Plates  4d,  4e).  However,  models  are 
prone  to  errors  in  the  strength  and  pathways  of  interior  in¬ 


ertial  jets,  larger  scale  SSH  biases  in  gyres  caused  by  errors 
in  subsurface  thermal  structure,  including  the  mixed  layer 
depth  and  temperature.  Such  errors  are  caused  by  the  atmos¬ 
pheric  forcing  as  well  as  by  the  model,  and  the  atmospheric 
forcing  is  a  significant  source  of  error  in  the  model  mean 
SSH,  for  example,  the  cleavage  in  the  large-scale  South 
Pacific  gyre  structure,  not  seen  in  either  observation-based 
mean,  is  caused  by  a  longstanding  problem  with  winds  from 
ECMWF.  These  winds  drive  a  South  Equatorial  Countercur¬ 
rent  that  is  too  strong  and  extends  too  far  east  [Metzger  et  al., 
1992],  a  problem  that  continues  in  the  1/12°  HYCOM  simu¬ 
lation  with  ECMWF  ERA- 1 5  forcing.  In  the  North  Atlantic, 
it  drives  a  western  boundary  current  that  is  too  weak. 

Because  of  the  preceding  problems,  it  is  necessary'  to  make 
corrections  to  the  model-based  means  using  ocean  data,  in¬ 
cluding  observation-based  means,  mean  frontal  analyses 
based  on  satellite  infrared  imagery,  SSH  variability,  the  limi¬ 
ted  number  of  nearly  simultaneous  air-dropped  expendable 
bathythermograph  (XBT)  under-flights  of  altimeter  tracks 
(very'  useful  for  the  Gulf  Stream;  Chassignet  and  Marshall 
[2008]  this  volume),  and  the  limited  number  of  often-repeated 
high-resolution  lines  of  XBT  or  acoustic  Doppler  current 
profiler  data.  This  was  done  to  a  limited  extent  in  the  HYCOM 
mean  used  in  the  data-assimilative  1/12°  global  HYCOM 
(but  not  the  mean  shown  in  Plate  4).  In  addition,  the  best 
available  atmospheric  forcing  is  chosen,  and  corrections  are 
made  to  the  wind  field,  for  example,  unrealistic  wind  curl 
features  around  Hawaii  were  corrected  to  eliminate  unrealis¬ 
tic  zonal  currents  w  est  of  Hawaii,  and  corrections  have  been 
made  to  avoid  having  wind  values  over  land  from  coarser 
resolution  atmospheric  fields  force  sea  points  in  a  finer  reso¬ 
lution  ocean  model  [Kara  et  al.,  2007].  Scatterorneter  data 
could  also  be  used  in  making  corrections  to  the  wind  speeds 
from  atmospheric  prediction  centers. 

2. 7.  Atmospheric  Forcing,  a  Significant  Data  Type  in  Ocean 
Prediction 

In  the  preceding  subsection,  atmospheric  forcing  play  ed  an 
essential  class  3  role  in  helping  eddy-resolving  1/12°  global 
HYCOM  simulate  an  accurate  mean  SSH  field  for  use  in  al¬ 
timeter  data  assimilation.  In  section  2.5,  it  was  shown  that 
forecasts  of  SSH  in  the  Persian  Gulf  and  at  coastal  tide  gauge 
locations  are  sensitive  to  the  quality  of  atmospheric  forc¬ 
ing  due  to  class  1  responses  in  shallow  water  and  in  coastal 
and  equatorial  wave  guides.  In  a  model  without  ocean  data 
assimilation,  Zamudio  et  al.  [2006]  showed  that  verifiable 
eddies  can  be  generated  by  a  combination  of  Tehuantepec 
winds  (class  1)  and  coastal  trapped  waves  originating  in  the 
equatorial  Pacific  (class  4),  including  subsequent  class  4  eddy 
propagation  -1 ,000  km  offshore.  In  another  non-assimilative 
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Plate  5.  Daily  observed  buoy  SST  (black  lines)  versus  1/12°  global  HYCOM  (red  lines)  using  three-hourly  FNMOC  NOGAPS 
wind  and  thermal  forcing.  The  time  series  span  from  August  1  to  September  30,  2005  and  cover  the  passage  of  Hurricanes  Katrina 
and  Rita  at  (a)  NDBC  buoy  42040  south  of  Mobile  Bay  and  (b)  buoy  42036  southeast  of  Pensacola.  FL. 
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Plate  6.  (a-f)  Verification  of  SST  and  altimeter  track  SSH  assimilation  in  1/12°  Gulf  of  Mexico  HYCOM  using  NCODA  and  data 
simulated  by  the  model  during  a  different  time  period  as  sampled  at  observed  SST  and  altimeter  traek  locations  during  that  time. 
Mean  error  (red)  and  rms  error  (blue)  versus  depth  to  500  m  for  (a,  d)  temperature,  (b,  e)  salinity,  and  (e,  f)  zonal  component  of  ve¬ 
locity  at  the  beginning  of  the  assimilation  period  on  August  29,  1999,  and  50  days  later  on  October  1 8,  1999.  (g-i)  Results  of  assimi¬ 
lating  real  SST,  temperature  profile,  and  altimeter  track  data  into  1/12°  Gulf  of  Mexico  HYCOM  and  verification  against  operational 
NAVOCEANO  frontal  analyses  of  satellite  infrared  imagery  overlaid  on  HYCOM  SSH  fields  for(g)  October  18,  1999,  (h)  Mareh  6. 
2000,  and  (i)  May  4,  2000.  Black  dotted  segments  of  the  frontal  analyses  arc  based  on  infrared  imagery  more  than  4  days  old. 
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model  simulation,  Zamudio  et  al.  [2002]  demonstrated  that 
Hurricane  Juliette  drove  coastal  downwelling  (class  1 ),  which 
then  propagated  along  the  coast  and  eyclonically  around  the 
Gulf  of  California  as  a  coastal  trapped  w/ave  (class  4).  Both 
of  these  are  examples  where  the  generation  was  a  class  1  re¬ 
sponse  to  atmospheric  forcing,  which  then  became  a  class  4 
feature  that  propagated  away  from  the  region  of  origin 

Simulation  of  SST  and  MLD  are  examples  where  a  combi¬ 
nation  of  class  1  and  class  3  responses  are  significant  [ Kara 
et  al .,  2003;  Kara  and  Hurl  hurt,  2006;  Kara  et  al..  An  eddy¬ 
resolving  oeean  model  for  the  Paeific  Ocean:  Part  2  :  Daily  and 
monthly  SST  variability  over  the  Pacific  Ocean  from  1 990  to 
2004,  submitted  t  o  Journal  of  Geophysical  Research ,  2007; 
hereinafter  referred  to  as  Kara  et  al.,  submitted  manuscript, 
2007],  class  1  for  the  sub-monthly  fluctuations,  and  class  3 
for  the  seasonal  eyele  and  its  interannual  variations  [Kara  el 
al .,  2004].  With  no  ocean  data  assimilation,  HYCOM  and 
NLOM  simulate  daily  time  series  of  SST  with  a  median  rms 
error  of  ~0.8°C  in  comparison  to  year-long  daily  time  se¬ 
ries  of  moored  buoy  data  [Kara  and  Hurl  burl,  2006;  Kara 
et  al.,  submitted  manuscript,  2007],  However,  with  assimila¬ 
tion  of  SST  from  advanced  very  high-resolution  radiometer 
(AVHRR)  data  in  the  form  of  model-independent  analyses, 
the  median  nowcast  SST  error  in  1/16°  and  1/32°  global 
NLOM  decreased  to  0.35°C  [Smedstad et  al .,  2003;  Shriver  el 
al .,  2007]  when  verified  against  unassimilatcd  daily  buoy  time 
series,  data  whieh  HYCOM  assimilates.  While  the  assimilation 
of  SST  and  SSH  data  affects  the  MLD  in  a  data-assimilative 
model,  their  assimilation  is  insufficient  to  accurately  constrain 
MLD,  and  the  quantity  of  subsurface  data  is  insufficient  to  per¬ 
form  that  role.  As  a  result,  there  is  a  heavy  burden  on  the  model 
to  use  a  combination  of  atmospheric  forcing  and  ocean  data 
assimilation  in  nowcasting  ML  D. 

In  Plate  5,  daily  SST  time  series  from  1/12°  global  HY¬ 
COM  (without  ocean  data  assimilation)  are  compared  with 
Gulf  of  Mexieo  buoy  data  from  two  locations  during  August- 
September  2005.  Both  show  realistic  model  response  to  the 
passage  of  Hurricanes  Katrina  and  Rita.  Global  HYCOM 
simulates  the  observed  drop  in  SST  after  the  passage  of  both 
storms,  indicating  realistic  upwelling  and  mixing  of  sub¬ 
surface  water,  as  well  as  sufficiently  accurate  atmospherie 
wind  and  heat  flux  forcing  to  achieve  this  result.  This  is  an 
important  capability  for  an  eddy-resolving  ocean  model  that 
could  make  it  very  useful  in  coupled  air-ocean  prediction  of 
hurricane  intensity. 

2.8.  Use  of  Simulated  Data  From  an  Eddy- Resolving  Ocean 
Model  in  Testing  Ocean  Dal  a  Assimilation 

Another  significant  use  of  eddy-resolving  ocean  models  is 
to  provide  realistie  simulated  data  for  a  variety  of  applica¬ 


tions,  as  discussed  in  section  1.3.  Providing  simulated  data 
for  testing  ocean  data  assimilation  techniques  is  a  very  use¬ 
ful  example  in  eddy-resolving  global  ocean  prediction.  In 
addition,  the  data-assimilative  results  can  be  used  in  observ¬ 
ing  system  simulation  and  assessment.  Here,  simulated  data 
from  1/12°  Gulf  of  Mexico  HYCOM  run  with  six-hourly 
interannual  forcing  was  used  in  testing  the  ability  of  the 
NCODA  system  [Cummings,  2005]  to  reconstruct  the  3D 
fields  of  temperature,  salinity,  and  velocity  when  only  SSH 
along  altimeter  tracks  and  SST  were  assimilated.  Complete 
model  fields  from  the  simulation  were  available  to  verify  the 
results  of  the  assimilation. 

The  3D  NCODA  MVOl  oeean  analysis  variables  include 
temperature,  salinity,  geopotential,  and  velocity,  which  are 
analyzed  simultaneously  using  a  HYCOM  forecast  from 
the  previous  assimilation  eycle  as  the  first  guess.  In  support 
of  HYCOM,  a  new  analysis  variable  was  added  to  NCODA 
that  corrects  the  model  isopycnal  layer  pressures  based  on 
differences  between  the  depths  of  density  surfaces  predicted 
by  the  model  and  those  from  observations.  The  NCODA 
horizontal  correlations  are  multivariate  in  geopotential  and 
velocity,  thereby  permitting  adjustments  to  the  mass  field  to 
correlate  with  adjustments  to  the  velocity  field.  The  velocity 
adjustments  are  in  geostrophic  balance  with  the  geopotential 
increments,  and  the  gcopotential  increments  arc  in  hydro¬ 
static  agreement  with  the  T  &  S  increments.  Synthetic  T  &  S 
profiles  obtained  using  the  approach  of  Cooper  and  Haines 
[1996]  were  used  for  downward  projection  of  SSH  below 
the  mixed  layer.  HYCOM  is  updated  by  adding  a  fraction 
of  the  analysis  increments  over  a  ehosen  number  of  time 
steps. 

In  the  test,  the  HYCOM  simulation  was  sampled  along  the 
observed  altimeter  tracks  of  TOPEX/Poseidon  and  ERS-2 
and  at  observed  locations  of  multi-ehannel  SST  data  during 
the  period  of  the  assimilation,  August  29,  1999  October  18, 
1999.  The  initial  condition  for  the  simulated  data  assimila¬ 
tion  experiment  was  the  simulated  state  a  year  later,  that  is, 
August  29,  2000,  when  the  Gulf  of  Mexieo  Loop  Current  and 
eddies  were  in  an  extremely  different  configuration,  for  exam¬ 
ple.  deep  penetration  of  the  Loop  Current  on  August  29,  1 999, 
compared  to  almost  no  penetration  on  August  29,  2000  Plates 
6a-6c  show  the  domain-averaged  rms  error  of  temperature, 
salinity,  and  the  zonal  component  of  velocity  versus  depth  at 
the  beginning  of  the  assimilation  period  (August  29,  1999) 
and  Plates  6d-6f  at  the  end  (October  18,  1999).  By  the  end  of 
the  assimilation  period,  the  error  has  been  greatly  reduced. 

The  test  was  then  repeated  with  assimilation  of  the  real 
altimeter  and  SST  data,  and  with  the  addition  of  temperature 
profile  data.  The  test  was  initialized  from  the  atmospherieally 
forced  simulation  on  August  29,  1999  and  run  for  1  year. 
SSH  analyses  from  the  HYCOM/NCODA  assimilation 
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were  compared  with  the  operational  NAVOCEANO  frontal 
analyses  based  on  AVHRR  imagery  (dashed  lines  overlaid 
on  the  SSH  fields;  Plates  6g-6i).  These  comparisons  indicate 
the  success  of  the  assimilation  in  mapping  mesoscale  vari¬ 
ability  in  the  Gulf  of  Mexico.  See  Chassignet  et  al.  [2005] 
for  a  comparison  of  five  different  data-assimilative  systems 
in  mapping  the  Gulf  of  Mexico  Loop  Current  and  eddies  in 
comparison  with  ocean  color  imagery,  including  an  early  ver¬ 
sion  of  1/12°  Atlantic  HYCOM  with  a  different  data  assimi¬ 
lation  approach,  the  1/16°  and  1/32°  global  NLOM  systems, 
1/8°  global  NCOM,  and  a  1/24°  Intra- Americas  Sea  NCOM 
In  other  testing,  HYCOM  has  been  used  to  improve  the  per¬ 
formance  of  error  covariances  in  NCODA,  and  ocean  models 
play  a  crucial  role  in  determining  error  covariances  in  some 
data  assimilation  methods,  such  as  ensemble  Kalman  filter¬ 
ing  [Evensen  and  van  Leeuwen ,  1996;  Brusdal  et  al.,  2003] 
and  ensemble  Ol  [Oke  et  al.,  2007]. 

3.  SUMMARY  AND  CONCLUSIONS 

Eddy-resolving  global  ocean  prediction  only  became  feasi¬ 
ble  near  the  turn  of  the  century  when  sufficient  data  and  com¬ 
puting  power  became  available  to  nowcast  and  forecast  the 
ocean  weather,  for  example,  the  surface  mixed  layer,  ocean 
eddies,  the  meandering  ocean  currents  and  fronts,  and  coast- 
ally  trapped  waves.  Computing  requirements  for  predicting 
the  global  ocean  weather  are  much  greater  than  for  atmos¬ 
pheric  weather  prediction  because  the  characteristic  diameter 
of  ocean  eddies  (-100  km)  is  20  to  30  times  smaller  than  for 
atmospheric  highs  and  lows.  Real-time  satellite  altimetry  is 
the  key  observ  ing  system  that  allows  global  prediction  of  the 
ocean  weather  associated  with  meandering  currents  and  ed¬ 
dies  (deep  ocean  mesoscale  variability),  but  sea  surface  tem¬ 
perature,  temperature  and  salinity  profiles,  and  atmospheric 
forcing  are  also  vital  real-time  data  sets  for  ocean  prediction 

Since  1998,  the  multinational  GODAE  has  fostered  the 
development  of  global  and  basin-scale  prediction  systems  of 
the  ocean  weather  in  several  countries  [Smith,  2000,  2006]. 
Because  of  limitations  in  computing  power,  most  of  these 
are  eddy-resolving  basin-scale  or  eddy-permitting  global 
prediction  systems.  As  of  2007,  only  two  eddy-resolving 
global  ocean  prediction  systems  have  been  developed. 

In  this  chapter,  we  have  focused  on  the  two  existing  eddy- 
resolving  global  ocean  prediction  systems,  especially  the 
critical  roles  that  ocean  models  play  in  these  systems,  and 
the  need  for  models  that  are  eddy-resolving,  not  just  eddy- 
permitting.  Resolving  the  first  internal  Rossby  radius  of  de¬ 
formation  has  been  used  as  a  criterion  to  distinguish  between 
eddy-permitting  and  eddy-resolving  models  because  of  the 
relation  of  the  Rossby  radius  to  the  predominant  space  scale 
for  baroclinic  instability.  This  criterion  is  useful  but  does  not 


provide  a  sufficient  definition  of  “eddy-resolving,”  and  ad¬ 
ditional  criteria  are  needed.  To  be  eddy-resolving,  a  model 
without  ocean  data  assimilation  needs  to  (1)  realistically 
simulate  the  eastward  penetration  of  inertial  jets  and  associ¬ 
ated  recirculation  gyres  where  they  exist,  (2)  demonstrate 
near-convergence  for  the  mean  strength  and  pathways  of 
both  upper  ocean  and  abyssal  currents  plus  their  mesoscale 
variability  statistics,  and  (3)  when  used  to  simulate  the  real 
ocean,  it  must  provide  realistic  comparisons  to  observations 
of  the  mean  upper  ocean  and  abyssal  circulation  and  their 
variability  statistics,  including  the  basic  physical  character¬ 
istics  and  behavior  of  eddies  and  meandering  currents.  Addi¬ 
tional  common  characteristics  of  eddy-resolving  models,  not 
seen  in  eddy-permitting  models,  are  a  basin-wide  explosion 
of  eddies  and  eddy-driven  abyssal  currents  constrained  by 
the  bottom  topography  which,  in  turn,  advect  the  pathways 
of  mid-latitude  upper  ocean  currents,  and  thus  influence  their 
mean  pathways.  This  upper  ocean-topographic  coupling  re¬ 
quires  that  mesoscale  variability  and  related  flow  instabili¬ 
ties  be  very  well  resolved  to  obtain  sufficient  coupling. 

In  an  ocean  prediction  system,  the  dynamical  interpolation 
skill  and  forecast  skill  of  an  eddy-resolving  ocean  model  are 
much  greater  than  in  an  eddy-permitting  model,  as  demon¬ 
strated  in  this  chapter.  In  most  regions,  an  ocean  model  with 
1/12°  equatorial  and  1/16°  (or~7  km)  mid-latitude  resolution 
for  each  variable  has  the  eddy-resolving  characteristics  de¬ 
scribed  above,  while  a  model  with  1/4°  or  1/8°  mid-latitude 
resolution  does  not.  A  doubling  of  the  mid-latitude  resolution 
to  1/32°  further  improves  the  dynamical  interpolation  and 
forecast  skill  of  the  model,  even  at  low  latitudes.  At  ~20°N, 
a  1/32°  model  with  assimilation  of  altimeter  data  from  three 
satellites  reliably  mapped  eddies  as  small  as  1/4°  in  diameter 
when  compared  to  an  unusually  clear  ocean  color  image. 

Model  dynamical  interpolation  skill  affects  the  quality 
of  the  nowcast  in  several  ways.  In  performing  a  nowcast,  a 
model  forecast  is  often  used  as  the  first  guess  for  assimila¬ 
tion  of  new  data.  Through  this  use  of  the  forecasts,  (I)  the 
model  helps  make  effective  use  of  delayed  data  and  helps 
fill  in  the  space-time  gaps  between  observations  (including 
temporal  gaps  between  the  repeat  cycles  of  satellite  altime¬ 
try)  by  using  model  predictive  skill;  (2)  the  model  helps 
convert  the  better  observed  surface  fields  into  subsurface 
structure;  (3)  it  helps  convert  the  better  observed  atmos¬ 
pheric  forcing  functions  into  useful  oceanic  responses,  and 
(4)  it  helps  apply  bottom  topography,  coastline  geometry, 
and  the  ocean  surface  as  physical  and  dynamical  constraints 
with  associated  boundary  layers.  As  additional  examples,  an 
eddy-resolving  model  should  also  play  a  major  role  (a)  in 
determining  the  mean  SSH  used  in  altimeter  data  assimi¬ 
lation,  (b)  in  separating  eddy-driven  steric  and  non-steric 
contributions  to  SSH,  and  (c)  in  testing  and  improving  data 
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assimilation  techniques  and  error  covariances  by  providing 
realistic  simulated  data. 

In  ocean  prediction,  it  is  essential  to  consider  the  classes 
of  oceanic  response  to  atmospheric  forcing  for  the  phenom¬ 
ena  and  regions  of  interest  (Table  1).  In  ocean  noweasting, 
this  is  important  in  considering  the  relative  impact  of  the  at¬ 
mospheric  forcing  and  different  oceanic  data  types.  In  ocean 
forecasting,  these  classes  impact  whether  or  not  the  time- 
scale  for  oceanic  predictive  skill  is  limited  by  the  timescale 
for  atmospheric  predictive  skill,  and  also  whether  a  coupled 
ocean-atmosphere  model  would  be  advantageous.  Results 
from  the  existing  eddy-resolving  global  ocean  prediction 
systems  demonstrate  forecast  skill  for  about  1  month,  glo¬ 
bally  and  over  most  subregions.  Outside  surface  boundary 
layers  and  shallow  water  regions  the  forecast  skill  typically 
is  only  modestly  impacted  by  reverting  toward  climatologi¬ 
cal  forcing  after  the  end  of  the  atmospheric  forecast  versus 
using  analysis-quality  forcing  during  the  entire  forecast.  The 
quality  of  the  atmospheric  forcing  during  the  forecasts  gen¬ 
erally  has  greater  than  average  impact  in  equatorial  regions. 

Surface  boundary  layers  and  shallow  water  regions  can 
respond  rapidly  to  the  atmospheric  forcing,  and  in  these 
locations,  the  timescale  for  oceanic  predictive  skill  can  be 
limited  by  the  timescale  for  atmospheric  predictive  skill.  In 
some  eases,  the  need  for  coupled  ocean-atmosphere  predic¬ 
tion  is  strongly  indicated,  for  example,  (1)  seasonal  to  inter¬ 
annual  prediction  of  phenomena  like  El  Nino,  (2)  hurricane 
prediction  where  SST  can  respond  strongly  to  the  hurricane 
and  the  hurricane  to  the  SST  change,  and  (3)  long-range  cli¬ 
mate  prediction  where  long  oceanic  adjustment  timescales 
can  be  important.  As  seasonal  to  interannual  and  long-range 
climate  predictions  become  increasingly  region-specific, 
eddy-resolving  ocean  models  will  be  needed  to  get  realistic 
current  pathways,  heat  and  mass  transports,  and  variability 
statistics,  all  of  which  are  essential  for  that  purpose. 

User  interest  in  real-time  ocean  products  has  been  high, 
as  evidenced  by  more  than  36  million  hits  to  the  NRL  Ocea¬ 
nography  Division  web  pages  during  2005.  Numerous  and 
diverse  research,  commercial,  military,  and  recreational  ap¬ 
plications  have  been  reported  [Hurlburt  et  a /.,  2002;  Sweet¬ 
s' tad  et  ai ,  2003;  Johannessen  et  a/.,  2006].  lee  prediction 
is  included  in  one  of  the  existing  eddy-resolving  global 
ocean  prediction  systems,  but  more  extensive  capabilities 
and  a  greater  range  of  applications  are  possible.  Increased 
capabilities  include  adding  tides,  more  advanced  data  as¬ 
similation  techniques,  ensemble  nowcasting  and  forecast¬ 
ing  (to  obtain  better  estimates  of  uncertainty  in  the  ability  of 
the  data  and  the  atmospheric  forcing  to  constrain  the  ocean 
model  nowcasts  and  forecasts,  and  to  help  identify  the  po¬ 
tential  for  extreme  events),  coupled  occan-atmosphere-icc- 
surfaee  wave  prediction  with  resolution  sufficient  for  hurri¬ 


canes,  global  biogcochemical  optical  tracer  and  contaminant 
prediction,  and  ultimately  earth  system  prediction  Many  of 
these  capabilities  would  greatly  increase  computational  re¬ 
quirements.  Based  on  the  results  presented  here  and  in  some 
of  the  references,  1/25°  equatorial  or  3-4  km  mid-latitude 
resolution  is  appropriate  for  a  pause  in  further  resolution  in¬ 
creases  for  global  ocean  prediction  systems  to  use  increases 
in  computing  power  to  add  such  capabilities,  and  tides  should 
be  included  in  a  1/25°  system.  Finer  resolution,  needed  in 
coastal  regions,  can  be  obtained  by  model  nesting. 
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